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ENTRENCHMENT OF THE WILLOW DRAINAGE 
DITCH, HARRISON COUNTY, IOWA* 


RAYMOND B. DANIELS 


Soil Conservation Service, U. S. Department of Agriculture, Ames, lowa 


ABSTRACI The changes in the Willow Drainage Ditch since 1919-1920 have been re 
constructed by comparison of the original profile of the ditch, historical records, and a 
survey of the ditch in 1958, Since construction, the drainage ditch has filled where the 
ditch was constructed on the Missouri River Valley. The drainage ditch in the Willow 
River Valley has entrenched, however, and becomes progressively deeper upstream, at 
taining a maximum depth of 42 feet at the Monona-Harrison County line. 

The drainage ditch in the Willow River Valley apparently deepened by channel 
scour and headward movement of knickpoints, but in a number of entrenchments rather 
than a single entrenchment followed by stabilization. Once a knickpoint has passed a 
point in the drainage ditch, stabilization of the channel] does not necessarily follow but 
channel scour may deepen the ditch more than passage of the knickpoint. 

The filling of the lower part of the drainage ditch probably is caused by the shary 
lecrease in gradient of the ditch, and by coinciding periods of high water in the Boyer 
Drainage Ditch or the Missouri River, to which the Willow is affluent. The entrenchment 
n the Willow River Valley probably was produced by a constructed increase in gradient 
of the drainage ditch as contrasted to that of the original river gradient in the same area 

The entrenchment of the Willow Drainage Ditch has been responsible for much of, 
but cannot explain all of, the entrenchment of its tributary streams. 


INTRODUCTION 


The rivers in Harrison County, lowa, commonly flooded their valleys dur- 
ing the early part of this century. To alleviate the flooding, the river systems 
were modified by construction of drainage ditches throughout the lower reaches 
of the rivers, Since construction, many of the drainage ditches have filled in 
the lower reaches, and at the same time have deepened and widened consider- 
ibly in the upstream reaches. 

The Willow Drainage Ditch, for example, has filled where the ditch was 
built in the Missouri River flood plain, but in the Willow River Valley it has 
deepened. Although the changes in the Willow Drainage Ditch are more 


severe than many of the other drainage ditches in Harrison County, these 


changes afford an excellent example of the method of stream entrenchment in 
the area. Also, the entrenchment of the Willow Drainage Ditch is important in 
explaining the entrenchment of its tributary streams. 

The Willow River heads in southwestern Crawford County, lowa, and 
flows in a southwesterly direc on through southeastern Monona County, In 
Monona County the Willow River is diverted into the Willow Drainage Ditch 
two miles north of the Monona-Harrison County line, The drainage ditch ex- 

Joint contribution of Soil Survey Investigations, Soil Conservation Service, U. S, Dept. 
of Agriculture; and the lowa Agricultural and Home Economics Experiment Station. 


Journal Paper No. J-3600 of the lowa Agrictural Experiment Station, Ames, Iowa. Project 
No. 1250. 


lo] 


Raymond B, Daniels—-Entrenchment of the 


R-44-W 
POTTAMATTAMIE COUNTY 


Location of Willow Drainage Ditch and the former channel of the Willow 
boundaries of the drainage districts do not necessarily coincide with the limits 
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tends from southern Monona County to southwestern Harrison County where 
it joins the Boyer Drainage Ditch south of Missouri Valley, lowa (fig, 1). The 
Willow Drainage Ditch is under control of the Harrison and Pottawattamie 
Drainage District and Upper Willow Drainage Districts No. 1 and 2 (fig. 1). 
The Willow River, upstream from the head of Upper Willow Drainage District 
No. 2, has not been straightened except in minor areas near bridges. 


CHARACTERISTICS OF THE WILLOW RIVER AND DRAINAGE DITCH 


Willow River Prior to Straightening 


Prior to straightening, the Willow River in Harrison and Monona Counties 
was a smal] stream that meandered considerably before flowing into the Mis- 
souri River Valley (fig. 1). The sinuous course of the Willow is shown when 
river miles are compared to valley miles, The stream length in Harrison County 
was approximately 26.3 miles, but the valley length is approximately 20.2 
miles. By comparison, the Boyer River in Harrison County, lowa, was approxi- 
mately 60 miles in length, and the Boyer River Valley 28 miles in length 
(Smith, 1888, p. 

Only sparse information is available on the width and depth of the river 
prior to 1900. Sales (1853, p. 244) recorded that in 1853 the channel of the 
Willow River in Monona County was 6 to 7 feet wide, 1.5 feet deep, and that 
the stream was flowing about 6 feet below the flood plain. 

In 1916-1918 the Willow River in Upper Willow Drainage Districts No. 1 
and 2 flowed 10 to 12 feet below the flood plain level, The bank-to-bank width 
of the stream ranged from a maximum of about 100 ieet to a minimum of 
about 60 feet.! No information could be found on the width and depth of the 
Willow River where it flowed into the Missouri River Valley. 

The gradient of the Willow River in Upper Willow Drainage District 
No. 1 averaged 5.2 feet per mile but ranged from a minimum of 1.7 feet per 
mile to a maximum of 7.9 feet per mile. In Upper Willow Drainage District 
No. 2 the gradient of the river averaged 7.5 feet per mile, and ranged from 
1.6 to 12.2 feet per mile (fig. 2). The wide range in gradients of the river in 
relatively short distances appears to be related to the entry of the larger tribu- 
taries into the river system, For example, the gradient of the channel upstream 
from Elk Creek (fig. 2, mile 15) is 6.4 feet per mile, but immediately below 
Elk Creek the gradient flattens to 1.7 feet per mile. Other examples of the 
flattening of the gradient below a tributary mouth are Thompson Creek, Huff- 
man Creek, and tributaries A, D, and E. The gradient of the Willow River 
downstream from the mouths of tributaries B, C, G, H, and South Willow 


Creek decreases slightly or is the same as the gradient upstream, An exception 


is the steepening of the channel downstream from tributary F (fig. 2). 

The flattening of the gradient of the Willow River immediately below the 
entry of a tributary into the river may have been due to some structural con- 
trol, but glacial till or bedrock are not known to be exposed in the channel of 
the drainage ditch. Generally, the sediments exposed in the drainage ditch to 

lhe width and depth of the Willow River were measured at the points that the Drainage 


Ditch crossed the river as shown by the plats and profiles of the drainage ditch on file at 
the Drainage Clerk's Office, Court House, Harrison County, Iowa. 
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depths considerably below the depths of the channel bottom of the original 
Willow River are dominantly silts with only sparse sand and gravel lenses. 
Thus structural control does not appear to have been responsible for the ma- 
jority of the inflection points in the Willow River. 

An alternative explanation for the flattening of the gradient below the 
tributaries is that the larger tributaries were overloading the Willow River 
with sediment and agerading the channel of the Willow for a short distance 
downstream. 

Price. in a survey of the Willow River in Drainage District No. 1 in 
1916. indicated that the river was filling (Harrison County, lowa, Drainage 
Record No. 7, p. 604°). Price suggested the filling was caused by the meander- 
ing of the river. 

The Willow River frequently flooded its valley to a depth of several feet 
is early as 1851, and the valley was, in general, considered unfit for cultiva- 
tion (Anderson, 185la. p. 190; 1951b, p. 168). In the early part of this cen- 
tury when the Willow River and its tributaries overflowed, they damaged crops 
on the inundated area. Local residents frequently disagree on the amount of 
the damage caused by flooding, but historical records indicate that the flooding 
frequently prevented planting or killed growing crops (Harrison County, lowa, 
Drainage Record No. 7, p. 604, 1104), The damage from flooding was espe- 
cially serious where the Willow River flowed in the Missouri River Valley 
(Harrison County, lowa. Drainage Record No. 3. p. 3). 


Characteristics of Constructed Drainage Channel 


lo alleviate the flooding by the Willow River, a drainage ditch was con- 
structed from Monona County, lowa, to the Boyer Drainage Ditch (fig. 1). 
Construction of the Willow Drainage Ditch in the Harrison and Pottawattamie 
Drainage District was started on June 20, 1906, and completed late in 1908 
(Harrison County, lowa, Drainage Record No. 3, p. 171, 198). 

The dimensions of the drainage ditch in the Harrison and Pottawattamie 
Drainage District were an 18-foot bottom, 1:1 side slopes, and a depth of 15 
feet from the top of the berm to the bottom of the channel®, The top of the 
berm was approximately 3 feet above the natural ground level. The total length 
of the drainage ditch was 6.6 miles. Total fall from the head to the mouth 
of the ditch was 13.5 feet, or an average of 2.04 feet per mile. 

The construction work on the drainage ditch in the Upper Willow Drain- 
ive District No. 1 started October 1, 1916. and was completed on October 16, 
1919. The ditch was cut with a bottom width of 12 feet, 1:1 side slopes, and 


15-foot berms. A strip of land 140 feet wide was condemned as right-of-wav 


for the drainage ditch. Average depth of cut was approximately 15 feet, but 
varied somewhat according to local conditions (fig. 3). The gradient of the 
ditch was 7.66 feet per mile throughout the drainage district (fig. 3), Total 
length of the ditch was slightly over 10.2 miles (Harrison County, Iowa, 
Drainage Record No, 7, p. 605, 683, 711). 


On file at the Drainage Clerk’s Office, Court House, Harrison County, Iowa. 


* The dimensions of the drainage ditch were obtained from the profile of the ditch on file 


n the Drainage Clerk's Office, Court House, Harrison County. Iowa. 
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Phe drainage ditch in the Upper Willow Drainage District No, 2 was con- 
structed between April 1, 1919, and June 1, 1920. The bottom width of the 
ditch was 12, 10, and 8 feet, with the width becoming less in the upper part of 
the district. Side slopes of 1:1 and a berm width of 15 feet were required in 
all segments of the drainage district, Depth of cut averaged 11 feet throughout 
the district (Harrison County, lowa, Drainage Record No. 7, p. 1106, 1140). 
The gradient of the ditch was variable, ranging from 12.14 to 8.45 feet per 
mile (fig. 3). 

In conjunction with construction of the drainage ditch in the Upper 
Willow Drainage Districts, most of the larger tributary streams also were 
straightened. An example of the amount of work on tributaries of the Willow 
is shown by the costs in the Upper Willow Drainage District No. 1. The esti- 
mated cost of the total project in drainage district no, 1 was $85,000, of which 
$49,630 was for construction of the main channel and the balance for the tribu- 
taries (Harrison County, lowa, Drainage Record No, 7, p. 689-690). 


Changes in the Willow Drainage Ditch Since Construction 


The Willow Drainage Ditch in the Harrison and Pottawattamie Drainage 
District (fig. 1) has been subject to filling since construction, At low water 
level the water flows 6 to 12 inches deep, and meanders over a bed of silt, The 
ditch was cleaned in 1916 and again in 1941, Prior to cleaning, the bottom of 
the ditch had filled with silt, and many of the drainage outlets were below the 
channel bottom (Harrison County, lowa, Drainage Record No. X, p, 469, 475, 
510). In 1954 Smith (Harrison County, lowa, Drainage Record No, XJ, p. 
834) reported that “ . . . Heavy sedimentation has filled up the bottom of the 
ditch greatly reducing its capacity. In places the bottom of the channel is only 
| or 2 feet lower than natural ground level.” Local residents believed that the 
Missouri River flood waters in April, 1952, had backed up the Boyer and 
Willow Drainage Ditches, causing considerable silting in the channel of the 
Willow Drainage Ditch and impeding drainage (Harrison County, Iowa, 
Drainage Record No, XI, p. 856). Inspection of the drainage ditch by the 
author in the fall of 1958 showed that the bottom of the ditch was 4 to 10 feet 
below natural ground level, and that numerous drainage outlets were partially 
or completely buried by the silt in the channel bottom. 

In direct contrast to the filling of the channel of the Willow Drainage 
Ditch in the Harrison and Pottawattamie Drainage District is the entrenchment 
of the drainage ditch in the Upper Willow Drainage Districts. 

Exact records on the changes of the Willow Drainage Ditch in the Upper 
Willow Drainage Districts prior to 1924 are not available, but interviews with 
local residents indicate that the ditch started to deepen and widen almost im- 
mediately, The earliest accurate record available is a measurement at a bridge‘ 
in drainage district no, 1 (fig. 3, mile 4.1) in 1924, The drainage ditch had 
deepened 6 feet at this point in approximately 5 years, Other records prior to 
1931 are not available, but a farmer (oral communication, Mr. Earl Thompson, 
Woodbine, lowa, 1957) living in the Willow Drainage Ditch area stated that 
by 1930 this ditch in the vicinity of Thompson Creek (fig. 1) had widened to 


* Bridge records were obtained from county and state highway official records, 
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a top width of about 60 to 70 feet. and had deepened to 20 to 25 feet below 
the flood plain level 

The Willow Drainage Ditch by 1931 had entrenched approximately 8 feet 
at mile 9.4, and had widened from an original width of 42 feet to about 57 
feet (fig. 3. cross section mile 9.4). The changes in the ditch at mile 5.6 were 
similar (fic. 3). In both cases the cross section of the channel was distinctly 
U-shaped, and not the trapezoidal shape of the original channel. 

County bridge records show that by 1936-37 the Willow Drainage Ditch 
had entrenched to about its present depth from the lower end of district no. 1 
upstream to mile 5.6 (fig. 3). Upstream from this point in the drainage ditch 
the records prior to 1940 or 1942 are not available because the bridges had 
been replaced. When the depth of the ditch in 1940 at mile 7.2 is compared 
to the depth at mile 8.3 in 1942, however, a distinct difference is noted (fig. 
3). Observations by Mr. Burkholder (oral communication, Mr. Verne Burk- 
holder, Logan, lowa, 1957) suggested that the ditch deepened slightly but 
widened considerably in the vicinity of mile 8.3 between 1941 and 1943, The 
sharp increase in gradient between miles 7.2 and 8.3 also suggests the drainage 
ditch was actively cutting in this segment between 1940 and 1942. 

The Soil Conservation Service, U, S. Department of Agriculture, con- 
ducted an Open Drainage Ditch Reconnaissance Survey of the Willow Drain- 
age Ditch in Harrison County in 1942. Records are available for parts of the 
drainage ditch. In the lower end of the survey the drainage ditch had deepened 
about 22 feet since 1920 (fig. 3, mile 12.5); upstream the ditch shallowed 
considerably and at mile 17.7 had increased the depth of its channel 11 feet 
since 1920 (fig. 3). 

lleff (1942, p. 31-41) had recorded a number of knickpoints on the tribu- 
taries of the Willow Drainage Ditch in 1942 from mile 12.5 upstream to mile 
17.7 (fig. 3). He also recorded that the bottom of the ditch was badly broken 
and checked in this segment, and many pronounced riffles with intervening 
stretches of slower water were present. Apparently the drainage ditch was ac- 
tively eroding between mile 12.5 and 17.7 just prior to 1942, This contrasts to 
the relatively stable channel bottom of the drainage ditch in most of drainage 
district no. 1 during the same period. 

(he channel of the Willow Drainage Ditch in 1958 is distinctly different 
than the original channel of the ditch in 1919-1920, The lower end of the 
drainage ditch in Upper Willow Drainage District No. 1 is 4 feet deeper and 


about 40 feet wider than the original channel, The ditch becomes progressively 


deeper and wider upstream, and has a maximum depth at the Monona-Harrison 
County line in Drainage District No. 2 (fie. 3, mile 19.6). At the Monona- 
Harrison Cou line the ditch has increased from an original depth of 11 feet 
in 1920 to a depth of 42 feet in 1958. The 1920 top width of 30 feet has in- 
creased to 110 to 120 feet 

Accurate records on the rate of deepening of the Willow Drainage Ditch 
ire not available. but a series of cross sections of the ditch at a bridge site on 
the north line of Sec. 16, T. 81N., R. 42W.. Harrison County, Lowa, illustrates 
the rapidity of entrenchment (fig. 3, mile 17.7). Between 1950 and 1952 the 


dit n h id deepened 15 reet ind widened about 20 feet. Since 1952 the bottom 
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PLATE 1 


\. Willow Drainage Ditch downstream from the mouth of Thompson Creek. The ditch 
this point is 36 feet deep and 110-120 feet wide. 


B. Knickpoint in the Willow Drainage Ditch, September, 1957, Note the action of 
water at the plunge pool. 


f the ditch has been relatively stable, but some adjustment in the steepness of 


the banks has taken place. 


The walls of the drainage ditch throughout the lower 17.7 miles of the 
ditch in the Upper Willow Drainage Districts are vertical, or nearly so, and 


the channel has a distinct U-shape (pl. 1A). Upstream from mile 17.7 the 
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1 top width of about 60 to 70 feet. and had deepened to 20 to 25 feet below 
the flood plain level 

The Willow Drainage Ditch by 1931 had entrenched approximately 8 feet 
at mile 9.4, and had widened from an original width of 42 feet to about 57 
feet (fig. 3, cross section mile 9.4). The changes in the ditch at mile 5.6 were 
similar (fig. 3). In both cases the cross section of the channel was distinctly 
U-shaped, and not the trapezoidal shape of the original channel. 

County bridge records show that by 1936-37 the Willow Drainage Ditch 
had entrenched to about its present depth from the lower end of district no. 1 
upstream to mile 5.6 (fig. 3). Upstream from this point in the drainage ditch 
the records prior to 1940 or 1942 are not available because the bridges had 
been replaced, When the depth of the ditch in 1940 at mile 7.2 is compared 


to the depth at mile 8.3 in 1942, however, a distinct difference is noted (fig. 
$). Observations by Mr. Burkholder (oral communication, Mr, Verne Burk- 
holder, Logan, lowa, 1957) suggested that the ditch deepened slightly but 
widened considerably in the vicinity of mile 8.3 between 1941 and 1943. The 


p increase in ¢ radic nt between miles 1.2 and 6.3 also suggests the drainage 


shar 
ditch was actively cutting in this segment between 1940 and 1942. 

The Soil Conservation Service, U, S. Department of Agriculture, con- 
ducted an Open Drainage Ditch Reconnaissance Survey of the Willow Drain- 
age Ditch in Harrison County in 1942. Records are available for parts of the 
drainage ditch. In the lower end of the survey the drainage ditch had deepened 
about 22 feet since 1920 (fig. 3, mile 12.5); upstream the ditch shallowed 
considerably and at mile 17.7 had increased the depth of its channel 11 feet 
since 1920 (fig. 3). 

Neff (1942, p. 31-41) had recorded a number of knickpoints on the tribu 
taries of the Willow Drainage Ditch in 1942 from mile 12.5 upstream to mile 
17.7 (fie. 3). He also recorded that the bottom of the ditch was badly broken 
and checked in this segment, and many pronounced rifles with intervening 
stretches of slower water were present, Apparently the drainage ditch was ac- 
tively eroding between mile 12.5 and 17.7 just prior to 1942, This contrasts to 
the relatively stable channel bottom of the drainage ditch in most of drainage 
district no. | during the same period. 

the channel of the Willow Drainage Ditch in 1958 is distinctly different 
than the original channel of the ditch in 1919-1920, The lower end of the 
drainage ditch in Upper Willow Drainage District No. 1 is 4 feet deeper and 


about 40 feet wider than the original channel, The ditch becomes progressively 


deeper and wider upstream, and has a maximum depth at the Monona-Harrison 
County line in Drainage District No. 2 (fig. 3, mile 19.6). At the Monona- 
Harrison County line the ditch has increased from an original depth of 11 feet 
1920 to a depth of 42 feet in 1958. The 1920 top width of 30 feet has in 
creased to 110 to 120 feet 
curate records on the rate of deepening of the Willow Drainage Ditch 
iailable. but ries of cross sections of the ditch at a bridge site on 
north line of Sec. 16, T. 8IN.. R. 42W.. Harrison County, Iowa. illustrates 
rapidity of entrenchment (fig. 3, mile 17.7). Between 1950 and 1952 the 


ch had deepened 15 feet and widened about 20 feet. Since 1952 the bottom 
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\. Willow Drainage Ditch downstream from the mouth of Thompson Creek. The ditch 
it this point is 36 feet deep and 110-120 feet wide. 


B. Knickpoint in the Willow Drainage Ditch, September, 1957, Note the action of 
water at the plunge pool. 


of the ditch has been relatively stable, but some adjustment in the steepness of 
the banks has taken place. 

The walls of the drainage ditch throughout the lower 17.7 miles of the 
ditch in the Upper Willow Drainage Districts are vertical, or nearly so, and 
the channel has a distinct U-shape (pl. 1A). Upstream from mile 17.7 the 
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ive L shape of the ditch is modified by abundant slump blocks, and the 
tions of the channel are very irregular. 


se 


The drainage ditch has been cut into a compact silt, and only sparse sands 

occul throughout il en h. The base of the stream. however. is flowing on silt. 

thick, that are constantly being shifted by the 

Above the Monona-Harrison County line the water is flowing directly 

1 compact silt, and only in the pools or quiet stretches is a thin silt, sand. 
ravel be 


The gra 


sand, and gravel, 1 t 


: drainage ditch in 1958 also differ from the gradients 
the upper Willow Drainage District No, 1 a uniform 
radient of 7.66 mile was constructed in the drainage ditch (fig, 3) 
The gradients in 1958 ranged from 4.7 to 6.7 feet per mile. In Upper Willow 
Drainage District No. 2 the original gradients of the drainage ditch ranged 
from 7.04 to 10.56 feet per mile, but in 1958 ranged from 5.1 to 21.2 feet per 
mile. The gradients of the ditch in 1958 are, in general, less than those of the 
original channel. but somewhat more than the original Willow River channel 
he, 2) 
ENTRENCHMENT 
Vethod 


Recorded Willow 


1! Drainage Ditch in the vicinity of the 
\onona-Harrison line since 1953 permit an analysis of the method of 

trenchment of the tch. 
\ knickpo nt 


e on the 


t in the Willow Drainage Ditch slightly above a 
9, T.8B1N., R.42W., Harrison County, Iowa, in 
The knickpoint moved upstream, and in November, 
ls north of the Monona-Harrison County line (oral 


unlap. owa, etween ovember. 
Dunlap, | 19 Bet N 


ward movement of the knickpoint was ap 
seen by the author on July 9, 1957. and had a 
Between July 9, 1957, and April 27, 1958, the 
ipstrea 5-90 feet and the vertical height increased to 
late 1B). From April 27 to May 1] 


600 teet: if 


int Move d 


, 1958, the knickpoint moved upstream 
was then essentially stationary until July 1, 1958, but be 
July 1 and August 15 moved upstream 1400 feet before disappearing 
ries ol riffl = 

The rapid movement of the knickpoint between April 27-May 1 and July 
Aucust 15, 1958 e related to periods of high water in the drainag 
tch because rainfall exceeding 1 inch occurred in the drainage basin during 
periods. The movement of the knickpoint could not be attributed to cut- 


1 more friable materials, such as sands, because the sediments were com 
t silts throughout the area of cutting. 


During the pel 


may D 


riods of slow headward movement in 1957 and 1958, the 
knickpoint was characterized by a vertical wall (pl. 1B). As the face of the 


knickpoint was undercut, shear planes developed in the unsupported sediments 
the l ( 


undercut sediments slumped, and a vertical face was reestablished, These 
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PLATE 2 


the plunge pool of a knickpoint. Note the wave action produced. 


B. Mass movement subsequent to passage of a knickpoint, Two distinct and one in- 
distinct slump blocks are visible in the upper right of the photo, The surface of the blocks 
represent former flood plain levels. 


processes were responsible for slow headward migration, The action of the 
water flowing over the knickpoint was similar to that described by Ireland, et 
al, (1939, p. 50, fig. 20). Some of the water always flowed down the face of 
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the knickpoint, even when undercut, but the percentage depended upon the 
stream discharge. A plunge pool 2 to 6 feet downstream from the base of the 
knickpoint always was present, and a swirling or wave action (pl, 2A) was 
produced in the plunge pool by the falling water. The wave action (pl, 2A) 
undoubtedly assisted in undercutting at all stages of the knickpoint develop- 
ment, Seepage water al the base of the knickpoint possibly was as important 
in undercutting as the back trickle and wave action of the plunge pool, The 
part of the knickpoint that undercut was usually at or below the seep line, and 
the saturated sediments possibly were less resistant to erosion than the drier 
overlying material. Thus, the water flowing down the face of the knickpoint, 
the swirling or wave action at the base, and the saturated sediments at the 
base probably contributed to the undercutting and slow upstream movement 
of the knickpoint. 

Once a knickpoint has moved upstream from a point in the drainage 
ditch, several adjustments in the channel take place. The passage of a knick- 
point normally exposes a seep line above the channel bottom, and mass move- 
ment of the saturated material is common (pl. 2B). The mass movement of 
the saturated material weakens the overlying sediments, and may be largely 
responsible for the formation of slump blocks. Slump blocks up to 30 feet in 
width have been measured downstream from an active knickpoint in the Wil- 
low Drainage Ditch, and as many as three distinctly separate slump blocks 
have been observed. The surface of each slump block was the surface of the 
original flood plain. In most cases the slump blocks have rotated backwards 
toward the valley walls. The blocks apparently move slowly downward toward 
the stream bottom and eventually are removed by the stream unless stabiliza- 
tion of the slumped mass occurs. 

Deepening of the Willow Drainage Ditch has not been entirely by head- 
ward movement of knickpoints, For example, between 1955 and 1956 the 
knickpoint in the drainage ditch moved upstream approximately one-third 
mile (fig, 4), but the channel continued to entrench downstream from the 
knickpoint. Little change occurred in the channel upstream, In June, 1957, 
heavy rains in the drainage basin resulted in a large flow of water in the ditch 
and deepened the channel considerably between mile 19 and 21 (fig, 4). Up- 
stream from mile 21] little change in the channel occurred. Between July, 1957, 


and August, 1958, the knickpoint continued to advance with concurrent deep- 
ening of the channel downstream, but little change in the channel of the drain- 
age ditch upstream from the knickpoint was measured (fig, 4). Thus, the 
deepening of the drainage ditch by the passage of a knickpoint is minor in com- 
parison to that produced by channel scour subsequent to passage of the knick- 
point. 


The changes in the Willow Drainage Ditch since 1936 show that entrench- 
ment of the drainage ditch essentially is a process of upstream migration of the 
eroding surface. The channel upstream from the eroding surface, while subject 
to some channel scour during periods of heavy flow, is relatively stable. The 
process of erosion of the Willow Drainage Ditch is not unlike that detailed by 
Ruhe (in press) for the process of pedimentation in southwestern Iowa, 


Raymond B. Daniels—Entrenchment of the 
Cause 


Since construction of the Willow Drainage Ditch, two changes have taken 

lace. The ditch has filled on the Missouri bottom, but has entrenched in the 

Willow River Valley. The gradient of the ditch in the Missouri River Valley is 
about 2 feet per mile, and immediately upstream in the Willow River Valley 
6.1 feet per mile (fig. 3). The abrupt change of gradient of the Willow Drain- 
age Ditch where it enters the Missouri River Valley probably results in a de- 
crease in velocity of the water with the subsequent deposition of part of the 
suspended sediment and aggradation of the channel. The deposition of sedi- 
ment 1s undoubtedly increased if the period of high water in the Willow 
Drainage Ditch coincides with the period of high water in either the Boyer 
Drainage Ditch or the Missouri River because the water backs up in the Wil- 
low Drainage Ditch under these conditions. It is doubtful if much sediment is 
deposited during low water because during these periods the ditch carries one 


foot or less of water. 


The Willow Drainage Ditch in the Willow River Valley was cut at a high- 
er gradient than that of the Willow River in the same segments. To illustrate, 
the average gradient of the drainage ditch was 7.66 and 8.48 feet per mile in 
Upper Willow Drainage Districts Nos. 1 and 2, respectively, The average gra- 
dient of the Willow River in the same segments was 5.18 and 7.50 feet per 
mile, The increase in the average gradient of the drainage ditch over the Wil- 
low River was about 1-2 feet per mile. 

lhe smooth, straight channel and increased gradient of the drainage ditch 
must have increased the velocity of the water, at least initially, and undoubtedly 


were responsible for the entrenchment of the ditch. 


E ffects 


The main reason for straightening the Willow River was to prevent 
damage to crops by flooding (Harrison County, lowa, Drainage Record No. 
7, p. 604, 1104). Interviews with local residents indicate the last flood on the 
Willow Drainage Ditch in the vicinity of Thompson Creek (fig. 1) was about 
1925 (oral communication, Mr. Earl Thompson, Woodbine, Iowa, 1957). At 
the upper end of the drainage ditch the last flood was in 1942 (oral communi- 


cation, Mr. M. Bendict, Dunlap, lowa, 1958). An example of the water-carry- 


ing capacity of the drainage ditch is shown by the absence of flooding in 
Upper Willow Drainage Districts during a period of heavy rainfall in June, 
1957. In a 24-hour period a series of rains equalling or exceeding 7 inches 
failed to fill the drainage ditch more than half full (oral communication, Mr. 
Earl Thompson, Woodbine, lowa, 1957). During the same rains, however, 
flooding did occur at or slightly below the lower end of Upper Willow Drain- 
age District No. ] 

The absence of flooding in most of the Upper Willow Drainage Districts 
for a number of years is interesting in light of a report filed by the commis- 


sioners and engineer that surveyed the Willow River in Upper Willow Drain- 
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age District No. 2 in 1918 (Harrison County, lowa, Drainage Record No, 7, 
p. 1106). “Judged from other straightened channels, for instance the Boyer 
[Drainage Ditch], this proposed straightened channel should wash out in a 
couple of years sufficiently to carry the maximum floods.” It is doubtful, how- 
ever, if the commissioners and engineer visualized a channel similar to the 
Willow Drainage Ditch in 1958. 


Influence on Tributaries 


Accompanying the entrenchment of the Willow Drainage Ditch has been 
the entrenchment of the tributaries of the ditch, Thompson Creek (fig. 1), for 
example, has entrenched at its mouth accordant with the channel of the Willow 
Ditch. The Trench of the creek deepens upstream for about 2 miles before 
shallowing where till is exposed in the channel bottom. In its headwaters 
Thompson Creek is cut to a depth below the adjacent flood plain of about 40 
feet. The entrenchment of Thompson Creek is typical of the larger tributaries 
of the Willow Drainage Ditch. The entrenchment in the mouths of tributaries 
such as Thompson Creek can be related directly to the entrenchment of the 
Willow Drainage Ditch. 

Shimek (1909, p. 298) described a number of creeks in Harrison County 
that had entrenched by 1909, which was 10 to 11 years before construction of 
the drainage ditch was completed in the Upper Willow Drainage Districts. 
Thus, not all of the entrenchment of the tributaries can be related to the en- 
trenchment of the Willow Drainage Ditch. 


SLMMARY 


Since construction of the Willow Drainage Ditch, the ditch has filled in 
the Missouri River Valley and simultaneously entrenched in the Willow River 


Valley. The method of entrenchment of the drainage ditch essentially has been 


1 process of headward encroachment of an erosional front on a channel that is 
relatively stable. 

rhe cause of filling of the drainage ditch in the Missouri River Valley 
probably is the sharp decrease in the gradient of the ditch in this segment 
ompared to the gradient of the ditch in the Willow River Valley, Concurrent 
flooding in the Boyer Drainage Ditch, or the Missouri River, causes water to 
back up in the Willow Drainage Ditch and also results in the deposition of 
sediment in the Willow Ditch. 

Che entrenchment of the drainage ditch in the Willow River Valley prob- 
ibly was caused by the increase in gradient of the constructed channel over 
that of the Willow River. 

Entrenchment of the Willow Drainage Ditch is responsible for much of 
the deep entrenchment of its tributary streams. 
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THE EFFECT OF SEDIMENT TYPE 
ON THE SHAPE AND STRATIFICATION OF 
SOME MODERN FLUVIAL DEPOSITS 


S. A. SCHUMM 


. 5. Geological Survey, Denver, Colorado 


ABSTRACT. The percent silt-clay in sediment appears to be indicative of some physical 
properties of that sediment. Silt-clay is defined as that material passing the 200 mesh sieve 
or smaller than 0.074 mm. The shape of a stream channel cross section is found to be in- 
versely related to the weighted mean percent silt-clay in channel banks and floor, The 
shape of paleochannels and their included channel-fill deposits may be similarly related 
to sediment type. 

Observations along actively aggrading streams suggest that the stratification of chan- 
nel-fill deposits may also be related to the percent silt-clay in these deposits, Stratification 
planes may be concave in channel fills composed of predominant silt-clay; whereas, they 
may be horizontal where the silt-clay content of the alluvium is small, and even convex in 
those channel fills with negligible silt-« lay. 


INTRODUCTION 
During the progress of investigations by the U. S. Geological Survey into 
the mechanics of erosion and deposition in semiarid valleys and into the rela- 


tionship between stream channel shape and sediment type, some observations 
were made which may be pertinent to the interpretation of fluvial deposits, 


either of recent deposition or in paleochannels. 

Students of fluvial morphology have recognized that the type of sediment 
transported by a stream exercises an important influence on the character of 
that stream. In this paper the importance of sediment type to channel cross 
section is demonstrated, and it is suggested that both the shape of a fluvial de- 
posit and the stratification developed during aggradation are dependent to a 
large extent on the physic al characteristics of the alluvium. 


METHODS OF INVESTIGATION 


Sediment samples of bank and channel materials were collected and 
measurements of stable channel cross sections were made at numerous stations 
long streams in Colorado, Wyoming, South Dakota, Nebraska, Kansas and 
New Mexico. A more detailed discussion of the methods and selection of sec- 
tions will be given in a report, dealing with the effect of sediment type on 
stream morphology, under preparation by the writer (Schumm, in press a). It 
is sufficient here to state that none of the channels contained more than 40 per- 
cent gravel. 

The shape of each channel is expressed as a width-depth ratio, The depth 
of the channel was measured from the edge of the first permanent surface or 
bank encountered above the channel floor to the lowest part of the channel 
floor. It is, therefore, a maximum depth. Width was measured from one edge 
of this bank to a corresponding elevation on the opposite side of the channel. 

Composite sediment samples were taken of the channel floor and both 
banks. Only the surface inch of the channel and banks were sampled along the 
selected cross sections, Depending on channel width, samples of channel sedi- 
ment were taken at 10 to 20 points across the stream. 
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In the laboratory the samples were subjected to a standard grain-size 
were first sieved, and if they contained more than 20 


inalysis. The samples 


percent silt-clay they were then prepared for hydrometer analysis. One im- 
portant deviation from the usual procedure was made to aid in the dispersion 
he samples (personal communication, I. S. McQueen and R. F. Miller). 
was placed in 700-800 ml of distilled water and 


of t 

Fifty grams of the sediment 

illowed to stand overnig| Che liquid was then decanted, removing most olf 

the salts that might prevent total dispersion of the sample by sodium hexameta- 

phosphate. Various data concerning texture and sorting characteristics wer 
ve grain-size curves of each sample. 

roblems to resolve in studies pertaining to sedi- 

arameter representative of the physical proper- 

physical characteristics of sediments of various 

n some detail by Bourcart (1941), usually the 

tion of sediments has tended to obscure the transitional 

boundaries hetween sediment classes and the physic al charac ter- 

istics of each sediment class However. work by engineers in the field of soil 

mechanics has revealed the 

the general properties of the whole, For example, Burmister (1952) demon- 


importance of the finer fractions of sediment to 


strates that the frost-heaving and drainage characteristics, as well as the capil- 
larity of a sediment ited to the amount of silt and clay in sedimentary 
material, In izen’s effective size (D,o), that size for which 10 per 
lated to the permeability of a sediment, 
he relation of stream channel shape to sediment 
m channel, neither median-grain size nor any 
size could be correlated with channel shape; 
bank and channel sediments was related to 
dth-depth ratio. Percent silt-clay is defined as 
er than 0.074 mm, or that material passing the 
of 0.074 as a critical value is probably misleading 
the third decimal place. It is possible that the 
tical zone between 0.05 and 1.0 mm would be 
to use that of the 200-mesh sieve, which is gen 
ictical limit of sieving. 
els sampled, the bank material differed appreci 
nels. However, by expressing percent silt-clay as a 
inks and channels, it is possible to obtain a 
of sediment forming the perimeter of a 


bank and coarse channel sediments or with 


stream channel, either with fine 
essentially similar bank and channel sediments. In order to obtain this weighted 
mean value, the percent silt-clay in the composite channel sample is multiplied 
by channel width. To this value is added the percent silt-clay in bank samples 
multiplied by twice the channel depth. This total is divided by the channel 
width and twice channel depth. In the remainder of this paper the weighted 
mean percent silt-clay will be designated as M and will be so used throughout 


ne text 
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THE SHAPE OF STREAM CHANNELS AND FLUVIAL DEPOSITS 
The shape of each cross section expressed as a dimensionless width-depth 
ratio is plotted against M in figure 1 for 69 stable cross sections. The correla- 
tion between channel shape and sediment is such that 
F 
where F is the channel shape expressed as a width-depth ratio and M is the 
weighted mean percent of silt and clay in the bank and channel samples. 
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Fig. 1. Relation between width-depth ratio (F) and weighted mean percent silt-clay 
(M). The regression line was determined graphically, Standard deviation is 0.487 log units, 
ind correlation coefficient is 0.91. 


The wide range of channel shapes comprising those plotted on figure 1 
suggest that for alluvial channels in which pebbles and gravel cover only 
negligible parts of the channel, the shape of the channel is dependent on sedi- 
ment type expressed as M. It is interesting to note that for those sections near 
caging stations, the mean annual discharge ranges from 20 to 4,920 cfs. In 
addition, many of the streams are ephemeral, in contrast to the perennial flow 
of the Smoky Hill, Kansas , and Republican Rivers, Drainage area ranges from 
about 1.7 square miles for the smallest ephemeral stream to 56,710 square miles 
above the Kansas River at Topeka. In spite of the wide range in stream regi- 
men and sediment type, the correlation coefficient for the regression line as 
obtained graphically is 0.91. 


++ + + + + + + + + ++ 
+ Me fe + + +—+++4+++ 
| 
+ + + 444 + + — 
ee 
| 
100.90 — 
+ +4++4 + — + 4++4++ 
+ + ++ - ++ 
> + + + ++ +++ + 
+ +—+ +++4+4 a @ 
> + + + ++ > _ 
| iil wee 
+ + + + ++ - + 4 44 
| | a 
100 | 
+ +4+4+ + 
} + + + + + 
+> + + + > > + + + + 
+ > + + + + + 
\ 


The Effect of Sediment 7 ype on the 
rreiatior » vindicate the selection of 
| sed ment ch iter 
hat cl 


M as a parametel 
In addition, the correlation of figure 1 shows 
little silt 
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rents in the same sandy material, however, it is “semicircular in section as a 
result of constant slumping of the walls under water” (McKee, 1957, p. 133). 

On the basis of these experiments, it seems possible to form submarine 
channels with a low width-depth ratio in sandy alluvium or alluvium with low 
M. In this event, the environmental factor would be dominant, Nevertheless, in 
View of the relationship of figure l. the cross-sectional shape of paleochannels 
and their included channel-fill sediment may be dependent to a large extent on 
the percent silt-clay in the banks and channel of the ancient stream. 


THE STRATIFICATION OF MODERN FLUVIAL DEPOSITS 


In addition to its effect on the shape of stream channels, silt-clay has been 
found to exercise an important effect on the mechanics of sediment deposition 
and the resulting stratification of channel-fill deposits (Schumm, in press b). 
Stratification varies from concave to convex as silt-clay in the alluvium de- 
creases, 

The flume experiments of McKee (1957) led him to conclude that “the 
nature of cross-stratification formed by the filling of channels largely depends 
upon the position of water level.’ He shows that in those cases where a channel 
is filled during a period in which water does not overflow the banks of the 
channel, stratification is horizontal, but with water levels high enough to over- 
flow and completely submerge the channel and flood plain “deposits develop as 
curving strata that are concave upward.” His studies on the Colorado River 
delta (McKee, 1939) revealed that many small tidal channels show concave 
stratification. 


Field observations of actively aggrading channels of ephemeral streams 


ind the relationship shown in figure 1, formed the basis for the preparation of 


figure 2. In figure 2 four generalized cross sections of filled channels typical of 
those studied are presented. The channels shown in figure 2 were not observed 
in cross section; nevertheless, detailed studies of the downstream changes in 
channel shape and sediment character, as the channels were filled, give a good 
indication of the probable appearance of the channels in section. 

Differences in stratification appear to be related to differences in the 
:mounts of silt-clay in the alluvium in the study areas. The following discussion 
f active channel aggradation as studied in four areas, in which channel deposi- 
tion is typical of that sketched in figure 2, may help to explain why such 
stratihcation occurs, 

In the Sage Creek channel (fig. 2 A) the percent of silt-clay at the stable 
ross sections averages 67, Deposition occurs initially around slump blocks that 
have fallen into the channel, for in areas of highly cohesive sediment the 
slumped material may not be immediately swept away because of its resistance 
to disintegration. As deposition continues, it is greatest on the channel sides 
ausing narrowing of the channel (pl. 1-1). When the channel is filled suffi- 
ciently to permit overbank flooding, the patterns of cross stratification are very 
similar to those of McKee’s flume experiments and may be due, in part, to 


wverbank flooding. 
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PLATE 1 


1. View upstream toward a sharp bend on Sage Creek, Deposits averaging 60 percent 
silt-clay have accumulated on both the inside and outside of the bend. 

2. View upstream on Sand Creek, above reach of heavy aggradation, Channel sediment 
contains 32 percent silt and clay. 

3. View downstream on Sand Creek where channel aggradation is marked, Sand is de 
posited on the channel floor but is covered by a veneer of silt and clay deposited during 


flood recession, Lateral deposition has formed inner banks along both sides of the channel. 


4. View downstream on Sand Creek about half a mile downstream from photograph 1-3. 
The sediment which nearly fills the channel contains 91 percent silt-« lay. 


In the Sand Creek area (fig. 2 B) the percent of silt-clay at the stable 
channel sections is 20, As deposition begins the sand is deposited in essentially 
horizontal layers on the floor of the channel (pl. 1-2), and as this coarser ma- 
terial is deposited, the percent silt-clay increases downstream, The concentra- 
tion of finer sediment in a downstream direction increases until the channel 
sediment contains well above 50 percent silt-clay, and lateral deposition is 
noted in this reach of the channel (pl. 1-3). The two types of deposition con- 
tinue contemporaneously, finer on the banks and coarser sediments on the 
floor of the channel, but in the later stages of filling all the deposited material 
is silt-clay (pl. 1-4). Horizontally stratified sands are restricted to the lower 
part of the fill and progressively toward the center near the top of the fill (fig. 


2 B). 
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In the Arroyo Calabasas channel silt-clay is negligible, less than 10 per- 
cent; the channel fills from bottom to top, (fig. 2 C) and stratification should 
be relatively horizontal. Convex surfaces of deposition have also been noted in 
very wide sandy channels (fig. 2 D). The Medano Creek, Colo., channel is 
600 feet wide, but water was present during the time of survey only for a width 
of 100 to 150 feet across this channel. The water was flowing on a part of the 
channel cross section lying at an altitude of several feet above the lowest part 
of the channel. In wide sandy channels or on flood plains, where the water 
absorbed by the fill is great, that part of the channel on which water is present 
may be raised by deposition. The next flood may shift the position of low flows 
to the lowest part of the channel. The resulting fill will appear in section as a 
series of overlapping convexities unless planed flat by high flood flows (fig. 2, 
D). McKee (1957) has duplicated similar convex stratification in a laboratory 
delta tank (see his fig. 4). 

The stratification of channel-fill deposits, therefore, may vary not only in 
relation to the environment of deposition but also with sediment character. The 
above discussions of aggradation in relation to the percent silt-clay in channel 
sediments may be amenable to laboratory investigation, for the sequence of 
deposition could be studied by using flume sediments that contain variable 
quantities of silt-clay. 

The generalized cross sections presented in figure 2, are supported to some 
extent by observations made by others. For example, summarizing the litera- 
ture on the subject, Shrock (1948, p. 223, fig. 179D) states that in small shale- 
filled channels, “The shale laminae sag slightly in the axial part and rise gently 
toward the sides of the channel—deformation due in part to compaction.” This 
deposit resembles that of the Sage Creek channel and channel A (fig. 2). 

Studies in estuaries of the French Coast reveal that the fine fractions of 
sediment are deposited on the banks of the channels; whereas, the sand is de- 
posited on the floor of the channel (Bourcart, 1947). In addition, Denson’s 
(19506, p. 701) studies of uranium deposits led him to attribute the rise of the 
contact between the Shinarump member of the Chinle formation and overlying 
shale of the Chinle, above the axis of a paleochannel, to a compaction of finer- 
grained sediments deposited along the flanks of the channel. This may be 
similar to deposition in the Sand Creek channel and channel B (fig. 2). 

Actually, the channel-fill deposits are only a part of the fluvial sedimentary 
complex. After the channel has been filled, sediment-laden water spreads over 
the flood plain. Sediment, formerly confined to the channel, now is spread over 
the valley floor. The overbank deposits may initially conform to the broadly 


concave shape of the valley floor. Maximum deposition, however, will generally 


be concentrated near the channel, which is assumed to be located near the cen- 
ter of the valley. The surveyed cross sections of valleys, in which the channel 
has been completely filled and where flood plain aggradation is important, 
ire broadly convex. 


CONCLUSIONS 


From previous work it was concluded that M seems to serve as a parameter 
representative of the resistance to erosion or general behavior of sediment in a 


Schumm 


stream channel, containing only small amounts of gravel. Also, as M increases 
the width-depth ratio decreases according to the following relationship: 

35 
In general, it may be assumed that the shape of channel-fill deposits will ap- 
proach this relationship. That is, in fine sediments the channel will be relatively 
deep and narrow; whereas, in sediment containing less M the channel will be 
relatively wide and shallow. 

The type of stratification exposed in a channel-fill deposit may also be 
dependent to a large extent on sediment character, Fine-grained, cohesive sedi- 
ments will tend to adhere to the channel banks, forming concave stratification. 
Noncohesive sediments will be deposited directly on the channel floor, causing 
horizontal stratification in the fill. Convex stratification can be formed by a 
stream, flowing in a wide sandy channel as it builds up its bed, which shifts 
laterally during floods. Convex stratification results also from deposition over 
the surface of a flood plain after filling of the channel. 

MeckKee’s (1957) flume work suggests that the environment of deposition 
is aiso important, and therefore, in any one locality it may be necessary to de- 
termine which of the two factors, environment or sediment type, is of most 
importance in developing the shape and stratification of fluvial deposits, 

How far observations of modern deposits and active aggradation along 
streams can be extrapolated to other areas of deposition or to the geologic past 
is unknown, but it is hoped that other channel-fill deposits can be interpreted 
successfully in the light of these 


l 


limited observations. 
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ACCRETION-GLEY AND THE GUMBOTIL DILEMMA 


JOHN C. FRYE,’ PAUL R. SHAFFER,’ H. B. WILLMAN,* 
and GEORGE E. EKBLAW' 


ABSTRACT. Inclusion of materials of different origins under the terms “gley” and 
“cumbotil” has led to misinterpretations and to confusion, It is suggested that in Pleisto- 
cene geology the term “gley” be used only for the gleyed material in a developed soil 
profile, and the term “accretion-gley” be used to designate deposited materials. It is further 
suggested that the term “gumbotil” be restricted in application to materials developed in 
itu. If such a restriction is not generally accepted, the term “gumbotil” should be 
ibandon d 


GENERALITIES 


In the field of research directed toward the stratigraphy, classification, and 
chronology of continental glacial deposits, few concepts have exerted as far- 
reaching influence as the gumbotil concept. The term “gumbotil” was first 
introduced by Kay in 1916, and discussed in detail by Kay and Pearce in 
1920. Since that time workers in the north-central region of the United States 
have embraced the concept so completely that it might be described as a major 
element of dogma among Pleistocene stratigraphers. 

Gumbotil was defined (Kay, 1916; Kay and Pearce, 1920, p. 89) as 

a gray to dark colored, thoroughly leached, nonlaminated, deoxidized clay, very 
sticky and breaking with starchlike fracture when wet, very hard and tenacious 
when dry, and which is chiefly the result of weathering of drift. The name is in- 
tended to suggest the nature of the material and its origin. 
This definition is both empirical and genetic. As will be explained, it is this 
dualism of the definition that produces the present dilemma. 

During the same period when Kay and Pearce were developing the gum- 
botil concept, the modern form of soil science was developing in the United 
States. Perhaps the lack of integration between the genetic concepts of gumbotil 
ind the rapidly developing data and concepts of pedology was fortuitous, For 


several decades glacial stratigraphers have commonly identified and named 


gupbotil on the basis of the empirical elements of its definition, but have in- 


terpreted its significance on the basis of the genetic elements of the definition. 
During this period there were only rare attempts to expand and develop 
the genetic aspects of the gumbotil idea (Leighton and MacClintock, 1930), 
but the evidence that demonstrates dynamically developing soil profiles was 
being discussed widely by soil scientists, Many of the terms and interpretations 
used in explaining the morphology and genesis of soil profiles were derived 
from pioneering work in Russia. Among these imported terms, “gley” is par- 
ticularly important to our discussion because the empirical definition of gum- 
botil falls within the range of materials included under the term “gley”, 
Some pedologists have completely rejected Kay’s genetic interpretation of 
rumbotil. This is clearly stated by Krusekopf (1948, p. 1] >) as follows: 

All of the characteristics of the gumbotil clearly indicate that it is a water 
leposit, or lacustrine in origin. Gumbotil is always a surface formation, and is 
clearly differentiated from the yellow till upon which it rests, It was deposited as 
a gumbotil and did not become so subsequent to its deposition, It was weathered 
and leached when deposited, 
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On the other hand, field relations at some places strongly indicate that other 


deposits classed as gumbotil occur as a part ol a developed profile of weather 
ing as suggested by Kay in his original definition. 


The conflicting interpretations of gumbotil have been aggravated by the 

» 

broad use, in a descriptive sense, of the terms “gley” and “gleyed material” in 


science. This multiple use has been defined in the Soil Survey Manual 


1] 
as tollows: 


n a few inches below the surface of hydro 
continue on down for many 
irise through the gradual filling of a wet basin, 
E idded lo at the surtace and be ing gleyed be 
thick mass of gleyed material, which may be 
sandy types, Obviously the upper part belongs 
part does not. This illustration does not extend to 
G horizon is clearly a part of the 


feet essentially 


solum and has 


larv wil he ( 
More recently 


the use of the term “gley” 


to include clays of different 
rin has been illustrated in 


1 summary by Senstius (1958, p. 362) in which 
states in part: 


the h nd bluish or steel-gray color of subaqueous soils, 


gley; ... In this example gley 
pond or swamp... On the left side 
issociated with a podsol on dry land. 


using the same term, without modifiers, 


represents the 


of composit 


for a 
swamp or pond and also for a clayey material 


art of a soil profile. If the underlying material 
the term gumbotil likewise might have been 
modifiers 


eyed material formed largely by accretion, 
listinguished from developed gley, 


and that 
distinguished from deposited “gumbo 
ind their relationships in the field, Min 
ta determinable only in the laboratory 


this discussion is restricted to those data 


methods 


ts, while applying the term “gumbotil” in an 


ed the possibilities of diverse origins, Depos ts 


n of gumbotil have formed by: (1) the a 


i lakes ol the 


t finest constituents of the debris 
low accretion of fine-textured sediments in 
l| plain, the materials being derived from 
htly higher parts of the till-plain surface and 
1ccumulation in marshy environments of fine- 


ermixed with sheet-wash accretion, (4) soil-profile 


moved by sheet-wash 


textured, eolian sediments I 


deve lopme nt in a situation ol hig h water-table where prot ess of at cretion and 
in situ development a1 


e both operative, and (5) processes of gleization in a 
yped soil profile where 


initial low permeability of the parent material and 


topographic to create a poor drainage, Only the last of these 
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five situations meets the genetic requirements of the original definition of 
gumbotil. To illustrate the divergent usages of the term “gumbotil” a few 
selected examples are described from Illinois, lowa, and Kansas. 


EXAMPLES 
Examples of many of the various materials that have been properly and 
improperly classed as gumbotil are readily observable in Illinois. As preface 
to the description of a few selected examples from Illinois it is appropriate to 
observe that the extensive exposures of the Illinoian drift plain in coal strip- 
mines in western Illinois reveal at only a few localities material that meets all 


of the requirements for classification as gumbotil, even though many of the 


strip pits occur in flat topography well removed from valleys and though the 
Sangamon soil is developed on clayey till that overlies Pennsylvanian shale. 

Material that is excluded from the empirical definition of gumbotil only 
by its moderate calcium carbonate content occurs extensively at the top of 
l'azewell age till in southwestern Iroquois County and areas farther north that 
were covered by early stages of glacial Lake Watseka. A typical sequence is 
exposed in cuts in the SW, SW14 SW sec. 29, T. 24.N., R. 14 W., 2% 
miles north of Clarence. One foot of black surface soil overlies three feet of 
blue-gray clay and silt, massive to blocky, with a few scattered pebbles and a 
few small nodules of calcium carbonate. An irregular zone of silty sand and 
gravel separates this upper clayey layer from another massive bed of clay, silt, 
and pebbles which is underlain by typical till. 

In Ogle County examples of accretion-gley (formerly identified as gum- 
botil) occur in the area mapped by Shaffer (1956) as Wisconsin Farmdale 
drift. Typical of these are a road-cut exposure in the NE44 NE sec, 2, T. 23 
N., R. 9 E. (Shaffer, 1956, p. 22), and exposures in building excavations and 
cuts in the SW14 NEW SE sec. 9, T. 23 N., R. 8 E. At the latter locality 
auger holes penetrated at least 12 feet of accretion-gley without reaching till 
below. The material is massive, noncalcareous, silty clay, Dark irregular zones 
of humus occur in the upper 11% feet, and the remainder of the deposit is dove- 
gray to blue-gray with local, irregular mottling of tan and brown throughout. 
In the lower half there are some silty zones but no pebbles, The locality is at 
the crest of a shallow valley wall and within a mile of the front of the Shelby- 
ville drift. 

The fact that this material is accretion-gley rather than a developed clay 
in a soil profile is clearly attested in the field by the character and occurrence 
of the humic streaking in the upper part, by its uniformity throughout a 12- 
foot thickness, by the suggestion of rudimentary bedding in the lower part, by 
the absence of siliceous pebbles in an area of generally pebbly till, by the 
distribution of the silt content, and by the complete lack of soil structure. 

In the area where Illinoian till is the surficial drift, examples of accretion- 
gley that was once considered gumbotil are common in the vicinity of Kewanee, 
Henry County. Such a deposit, below thin Peorian loess, was observed in road 
eut and auger hole in the NE cor, SE sec. 26, T. 15 N., R. 5 E. The upper 


wo teet 


massive, 


ult 
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y. irregularly streaked with humus, that grades downward 
noncaleareous, silty clay, The underlying two 


ay. silty clay is distinguished by the presence of s« attered 


bbles. Below this the auger penetrated four feet of brown, pebbly, silty clay 


reous till. That at least the uppermost three feet is an ac- 


ittested in the field by the scarcity of pebbles, by the dis 
of the humus streaks, and by the texture. 

if the sharp localization of accretion-gley is exposed 

Vy NEY NW) sec. 36, T. 8 N., R. 2 E., Fulton 

downward succession consists of 10 feet of Peorian loess. 


rmdale loess, Sangamon soil, and IIlinoian till. A lens of 
1 maximum thickness of three feet occurs in a shallow de 
vian drift plain, This lens is sharply limited; it consists o! 
silty clay, with irregular humus streaks in the upper part. 
hough they occur in the adjacent Sangamon soil B-horizon. 
traced laterally under the accretion-gley lens, but it is 


s mottled with LV. suggesting some gleving of the weathered till 
le 


n deposit. Although the material in this lens meets the re 


definition of gumbotil, its depositional origin is 
lat onships. 

s lens the developed profile is similar to the Sanga 

less than a mile away and of which the min 


Brophy (1959), It contains a red-brown, clayey 


1 «le velo] d ciay skins, manganese pellets, and well preserved 


rades downward into tan-brown, calcareous Illinoian till. 


ie-cray, calcareous till with irregular zones of sand 


tains excellent examples both of accretion 
botil and of developed Sangamon soil profiles 
or gumbotil. Typical of the former is a road 
NEW sec. 5, T. 14 R. 4 W. Eight feet 
some tan mottling, massive, noncalcareous, 
the upper 14 feet is exposed below tan 
nches of weakly calcareous, tan silt, below 

feet of tan and gray. weakly calcareous 


of brown and gray, noncaleareous, silty 


side of a valley in the NW! { \W! { SW! j 
ses Sangamon soil, developed on Illinoian till. 
npirical and genetic requirements for gumbotil 


loesses the A-horizon, approximately two feet 


bles, and gray-tan in color, rests on six feet of 
which can he classed as sumbotil. The uppel 
s very clayey. has well developed soil structure, contains 
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mottling. The color grades downward through greenish gray to yellow-brown; 
the clay content, soil structure, and mangaense pellets and streaking diminish, 
and the sand content increases downward. Below the B-horizon is four feet of 
oxidized and leached till, five feet of oxidized and calcareous till, and four feet 
of blue-gray calcareous till. 

In Iowa, an example of accretion-gley on Kansan till is well exposed in a 
fresh road cut in the SEY, SW sec. 7, T. 82 N., R. 4 E., Clinton County. Be- 
low thin, tan loess is exposed 41% feet of gray, noncalcareous, silty clay resting 
on six feet of light gray, weakly calcareous till, grading downward to tan and 
vray, calcareous, Kansan till. The accretion-gley is irregularly streaked with 
humus in the upper part and has tan streaks in the lower part. Indistinct bed- 
ding with somewhat more silty zones is present in the lower part. The sleying 
effect of the local environment is illustrated by the blue-gray streaks and 
pendants extending as much as five feet into the till. 

Two localities in Kansas exemplify two types of gleyed material. A soil 
developed in Nebraska till and overlain by Kansas till is exposed in the NE%4 
SE sec. 6, T. 2 S., R. 20 E., Doniphan County, Kansas (Frye and Leonard, 
1952, p. 55 and 69). The soil has been classed as a gumbotil and is largely 
developed in situ from till parent material. It is leached, highly clayey, dark 
vray, massive, and in general will meet the requirements of the original defini- 
tion of gumbotil. Siliceous pebbles are present throughout although they are 
less abundant and smaller in the upper part, the contact of the B-horizon with 
the underlying material is gradational, the leached zone is little more than 
three feet thick, and the unleached till below the gley contains not only siliceous 
pebbles but also pebbles of limestone and igneous rocks. 

An example of accretion-gley occurs in upland situations in Atchison, 
Jackson, Jefferson, and Shawnee Counties in northeastern Kansas, In the NE'4 
sec, 12, T. 7 S., R. 18 E., Atchison County, a gleyed layer attains a thickness 

0 feet and has been informally called the “Nortonville clay” (Frye and 
onard, 1952, p. 81). This massive clay layer, which has been called gumbo- 
. is gray, free of pebbles and cobbles, and noncaleareous; it rests with a sharp 
contact on oxidized and cobbly Kansas till, and where not overlain by thin 
loess it contains a heavy planosol in its top. Its regional and topographic rela- 
tionships, its lack of siliceous pebbles and cobbles so common in the till im- 
mediately below, its uniformity throughout such a thickness, its sharp basal 
contact, as well as its chemical and mineralogical composition, clearly indicate 
that it is not the product of in situ weathering but is an accretion deposit on 
the initial surface of the Kansan till plain. 


CONCLUSIONS 


The diverse origins of deposits classed as gumbotil must be recognized, 


is has been done in the case of gley or gleyed materials, Our terminology 


needs clarification so that a deposit named empirically on the basis of physical 
features will not be automatically and unquestioningly given genetic signifi- 


cance, 
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We suggest that when used by Pleistocene geologists the term “gley” be 
restricted to clayey gleyed material in a developed profile and the term “ac- 
cretion-gley” be used for reference to the several types of slowly accumulating 
deposits of surficial clay. It is apparent that to be meaningful, the term “gum- 
botil’”” must be restricted to material resulting from in situ development, be- 
cause the indiscriminate use of this term for materials of different origins has 
produced more confusion than clarification, If such a restriction is not gen- 
erally accepted, it is our judgment that progress will best be served by abandon- 
ing the term “gumbotil.” Such deposits should then be described lithologically 
the same as other parts of the succession of sediments, It is appropriate in 


seological literature to use pedologic terms in describing a soil profile. 
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Rb-Sr ANALYSES AND AGE DETERMINATIONS 
OF CERTAIN LEPIDOLITES, INCLUDING 
AN INTERNATIONAL INTERLABORATORY 

COMPARISON SUITE 


LEONARD F, HERZOG*, WILLIAM H. PINSON, Jr.’ 
and PATRICK M. HURLEY 


ABSTRACT, Rb, Sr and radiogenic “Sr concentrations, and radioactivity “ages” have 
determined by stable isotope dilution for lepidolite samples from eight pegmatite 
currences 

Five of the samples analyzed comprise an “Intercomparison Suite” ranging in age 

om 100 < 10° to 2600 x 10° years. Portions of these samples were sent by us to several 
otl laboratories for analysis, as a test of the interlaboratory reproducibility of stable 
isotope dilution analysis. In this paper analytical data for this suite from four laboratories 
ire compared. Interlaboratory reproducibility and intralaboratory reproducibility are es- 
sentially identical for the age determinations; it is apparent that age ratios for lepidolite 
an be measured routinely with analytical accuracy better than 4 percent (standard 
leviation) 

Rb was determined with a precision of 1 to 2 percent, and Sr with a precision of 1 to 
3 percent within individual laboratories; difference between results obtained at the several 
laboratories may be due in part to small differences in Rb and Sr standard solution cali- 
brations and in part to sample inhomogenieties; the good interlaboratory agreement in age 
leterminations favors the latter interpretation. 

A comparison of stable isotope dilution, optical spectrographic, and neutron activation 
Rb/Sr data for sets of independently collected and analyzed samples from six localities is 
ilso given. For each pegmatite, there is good agreement (a range of a few percent) in the 
ages determined by isotope dilution and neutron activation, although Rb concentrations 
have ranges of up to 30 percent. 


The study gives insight into the sources of error in the optical spectrographic tech- 


niqt 


ie and provides standards which may make possible the improvement of this method of 


INTRODUCTION 

Phe Rb-Sr analyses of certain pegmatite lepidolites made by L. T. Aldrich, 

(Davis and Aldrich, 1953) by the stable isotope dilution method were in 
marked disagreement with earlier optical spectrographic analyses of samples 
from the same localities (Ahrens, 1949). To determine the source of the dis- 
igreement, at Prof. L. H. Ahrens’ suggestion our laboratory prepared an “In- 
tercomparison Suite” of five lepidolites (Herzog, et al., 1955). Portions of 
these samples were sent for analysis to seven mass spectrometer-equipped 
laboratories throughout the world which had expressed an interest in Rb-Sr 
dating. 

Dr. Aldrich and his colleagues at the Department of Terrestrial Mag- 
netism, Carnegie Institution of Washington have reported analyses of the five 
samples sent to them (Aldrich, Wetherill and Davis, 1956), E. Schumacher 

then at the Institute of Nuclear Studies of the University of Chicago) has 
reported an analysis of one sample in a letter to one of the authors (Schu- 
macher, personal communication, 1956), and Webster, Morgan and Smales 
1957) have published determinations of two samples made at Harwell, Our 
group has also analyzed all the samples in this suite at M.1.T. with the results 
Department of Geophysics and Geochemistry, The Pennsylvania State University, Uni- 
versity Park, Pennsylvania. 
* Department of Geology and Geophysics, Massachusetts Institute of Technology, Cam- 
bridge, Massachusetts 
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listed in t > 1 (Herzog. et 


il.. 1955. 1956: Hurley, 1956b). As far as we 


are aware, none of these samples has as yet been analyzed in any other labora- 


We have withheld publication of an analysis of this data for two years 
while aw iiling more analyses But the probability of other analyses of these 
samples being made in the near future now seems small, so we have decided 


to present the data at hand at this time. 


Des: ription of Samples 


The samples analyzed in this study ranged in grain size from fine powders 
of approximately 0.lmm grain diameter to well-crystallized books of up to 
uses, the book thicknesses were considerably (ap- 

less than their diameters. 


15mm diameter. In most 
pi ximately <4 to 10) 
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up of lepidolites with as wide a range of ages as 
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interested in o 
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ime that vas important 
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and geologic im 
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It is believed that all reporting laboratories analyzed the samples “as re- 
ceived’, i.e., with no additional treatment. 


ANALYTICAL TECHNIQUES 


The technique of isotope dilution has been described by several authors 
(Inghram, 1954; Hintenberger, 1956; Herzog and Pinson, in preparation). 
However, the analytical techniques used at the four reporting laboratories dif- 
fered in detail, and as these differences may be important in connection with 
differences in results, they have been summarized below. First, however, a 
summary description of our own procedures is given: 


M1.T. Nuclear Geophysic Section Procedure 

Lepidolite samples of 0.1 to 1.0 gram weight (table 1) were weighed out and digested 
with HF and HCI0,. Small amounts (up to 0.5 percent) of insoluble residue always re- 
mained after digestion. The residues from some samples were examined under the petro- 
graphic microscope and found to consist mostly of minute grains of topaz, Analysis by 
optical spectrograph indicated that they had no measurable Rb or Sr content. 

After digestion, each sample was mixed with a known volume of a calibrated “Sr- 
enriched strontium solution.’ Then it was diluted with HCl to a known volume (generally 
100 or 250 ml). 

In age determination work it is desirable to add spike to parent and daughter element 
simultaneously. Thereby one insures, as well as is possible, that any inadvertent losses of 
solution will affect each element proportionately to its concentration. In making these 
inalyses, however, we added “Rb and “Sr spike separately, adding Rb spike to only a 
fraction of the total sample, after it was in solution. This technique was adopted because 
of (1) the high ratio of Rb to Sr in lepidolite and our desire to (2) conserve spike solu- 
ons but (3) still use a large, representative sample. Lepidolite samples containing 
amounts of Sr appropriate for analysis contain more than an order of magnitude more Rb 
than can be accurately determined by the use of the maximum amount of Rb spike we felt 

was desirable to sacrifice for a single analysis.‘ Hence, an aliquot containing 1 to 5 
percent of the total solution was removed by volumetric pipette and placed in a “Vycor” 
sh, and a separately measured volume of concentration-calibrated “Rb spike solution 
vas added to it. The relative concentration of Rb in the aliquot was increased sufficiently 
for mass spectrometric analysis by the precipitation of Al (and Fe, if present) with am- 

m hydroxide; Sr was separated from the main portion of the sample by ion ex- 

techniques. Operation of the columns was monitored by optical spectro- 


j 
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Calibration of Sp ke Solutions The Rb and Sr spike solutions used in isotope dilu- 
analysis in this work were calibrated against standard chloride solutions of Rb and 
“normal” isotope composition, using the stable isotope dilution technique. Concen- 
ms of spikes were also determined by weighing them as they were made up, but as 
ike quantities weighed were small (less than a few milligrams, generally) it was 


t 


better accuracy could be obtained through the use of primary standards made 
ip with amounts of Rb or Sr salts large enough to be weighed very accurately. 
Large volumes of the primary standard normal-isotope-abundance Sr and Rb* solutions 
made up, so that they could be used for spike calibrations for many years.. However, 
solutions of normal Sr and Rb were also made up at various times as checks on the 


1 materials, in which the ratios of the various isotopes have been markedly alte red 


to those found in the naturally-occurring element, are commonly referred to as 
spikes” 


‘ The use of separate portions of the mineral sample for the analysis for each element is 
also possible. This procedure was rejected on the ground that it is most difficult to insure 
that both samples have the same composition, The differences between our results and 
those of other investigators in the Sr and Rb analyses listed below for the Intercalibration 


Suite, for some of the coarse-grained samples (e.g., Varutrask), indicate that such a pro- 


cedure might indeed have introduced appreciable errors in the determined Sr/Rb ratios, 
and, hence, in the determined ages. 
Because Rb is a congener of K, it is very difficult to obtain a pure Rb salt to use for 
the preparation of a standard solution, We used spectrographically-pure RbCl and purified 
further by ion-exchange resin techniques, testing for purity by spectrographic analysis. 
In a later experiment, one batch of “C.P.” RbCl was shown to contain 14 percent K, 
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ind Wetherill, 1956). Aldrich, et al. used about 50 mg of lepidolite per analysis, They 
chose to add Rb and Sr spike simultaneously, for this purpose using 150 to 300 micrograms 
of “Rb spike and 5 or 50 wgm of “Sr spike. 

Rb was placed on outgassed tantalum filaments as perchlorate and Sr as nitrate. The 
spectrometer used (built at D.T.M.) is nearly identical to the one used at M.LT, However, 
the isotope ratio measurement procedure used for Sr analysis was somewhat different: ion 
current was amplified by an Il-stage Allen-type electron multiplier (Allen, 1947) before 
being inserted into the Vibrating Reed Electrometer. To calibrate results obtained by using 
the multiplier, some samples were analyzed both “with” and “without” it. It was found 
that a correction of about 1 percent per mass unit to analyses made with the multiplier 

ld reconcile results in the Sr mass range. Aldrich, et al. believe that this effect is pri- 

y due to variations of the magnetic field in the multiplier region when the analyzer 
magnetic field is varied in peak “sweeping”. 


wou 


maril 


These investigators report a routine-analysis reproducibility for Rb isotope abundance 
inalysis of + | percent, and for Sr (because of isotope fractionation), up to 3 percent. 
Institute for Nuclear Studie {rgonne National Laboratory Procedure 

Dr. Schumacher, (Hurley, 1956s) used one of the 12” radius, 60° deflection, direction- 

sing spectrometers with electron multiplier at the Argonne Laboratory for his analyses. 

trument had a mass resolution® of 1 at m/e 800, compared to 150-250 for the 
used by the other analysts. An oxidized Ta filament was used, Rb was 

ide and Sr as oxalate. Some 50-60 peak sets were recorded per analysis. 
for mass discrimination was made in the determination of radiogenic 
on a measurement of the ratio “Rb/“Rb in normal Rb with the multi- 
m4 0.010 solution was made for a correction factor which would 
; s number agree with the accepted value, 2.590 + 0.003", obtained by Nier and 

ny Aldrich ind Herzog 

ay Sol 


d 


umacher used less than 0.5 micrograms of Sr ver spectrometer run, and 
I 


s than 0.05 nem of Rb. Pure spike experiments showed contamination levels of 0.1 
m Sr and 0.05 nem Rb, per analysis, for his entire procedure. 

ul Chemistry Group, A.E.R.E.) Procedure 

t al. (1957) report that their procedure is generally similar to that used 


Howe ver, the \ k R.E . worke rs added Rb spike to only a portion of the sample 


as was also done at M.I.T. Ion-exchange separation was done on Zeocarb-225 


id of the Dowex-50 used by other analysts. A 6” radius, 60° deflection, single- 
trometer was used: the detection circuit consisted of a type 954 
ind a D.C, amplifier, 


RESULTS 


M.I.T. analytical data for the Intercomparison Suite are presented in 

1. All the data necessary to recompute Rb and Sr concentrations and 

“aves” are given in this table. Our Rb, Sr and age determinations for the suite 

ire compared with those made by the other analysts, using portions of the same 

samples, in table 2. 

Rb concentrations in these lepidolites range from 0.56 to 3.3 percent by 

28.32 percent by weight of this Rb is *’Rb, the radioactive “parent” 

of “Sr. Sr concentrations vary from 6.5 to above 330 parts per million by 
weight, and nearly all of this radiogenic **Sr. 

Agreement in the *’Sr/**Rb ratio and apparent age columns is striking. 
All results agree within the standard deviations stated by the investigators, 
which are ordinarily in the range 3.5 percent to 7 percent. 

Davis, et al. (1953) had previously analyzed specimens collected inde- 
pendently from some of these localities. Revised data on these samples, given 
in recent papers (Hurley, 1956; Aldrich, 1956; Davis and Aldrich, 1956), are 
included in table 2. We have also reported in this table M.L.T. analyses of cer- 
tain lepidolites from localities which are not a part of the Intercomparison 


e., with the particular slit combination chosen. 


These values are erroneously printed as 8.583 and 8.590 in Hurley, 1956s. 
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Suite and, for comparison, analyses by Aldrich, et al. of other specimens from 


two of these localities. 

As a matter of interest Ahrens and Gorfinkle’s (195la) optical spectro- 
graphic analyses and Harwell neutron activation Rb analyses of samples from 
some of these localities (Webster, et al.. 1957) are also included in the table. 

The ages given in table 2 have all been computed using the same decay 
constant for **Rb, i.e.. A 1.385 X 107" disintegrations per year (cor- 
responding to a half-life of 50 X 10° years). This procedure of course does 
not affect the relative ages of the samples; the “ages” corresponding to any 
other half-life can be readily computed from the radiogenic **Sr/parent **Rb 


ratios given for each sample. 
Taste 1B 
Other Data Used in Computing Rb and Sr Concentrations and Ages 


fope abundance 
mass 84 § 85 nass $ mass 87 mass 88 
),0056 0.0702 0.9256 
0.4600 0.0890 0.3067 
0.7215 0.2785 


0.1017 0.2989 


pikes used (in micrograms per 
5. used for N45, Al07 
used for All0 
used for A106 and Al09 
61.0, used for N45 and Al07 
6.6, used for A106, A109. and All0 


c weights used in calculations 
tb/spike Rb 85.56,/86.80 0.9857 


63 /85.88 1.022 


.00/85.88 1.014 


spike St 87. 
Sy spike Sr 8 


‘ 


10° disintegrations, 


DISCUSSION OF RESULTS 
\ single lepidolite age determination by stable isotope dilution requires 
two analyses, one for Rb and one for Sr, and a standard solution and a spike 
solution must be prepared and calibrated for each element investigated, In the 
present study, each laboratory independently prepared the necessary standard 
and spike solutions, Also, at least in our laboratory, it was necessary to use 
more than one spike solution per element during the course of the investigation, 


because the original solution was used up. 


Intercomparison Suite 
{ges—In view of the number of operations involved in each age deter- 
mination, the very small ranges in the ages reported for particular samples 
ire evidence that isotope dilution is an analytical technique of high accuracy. 
L-N45.—The Bikita sample was analyzed by all four reporting labora- 
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TABLE 
Rb and Sr concentration, and age determinations reported from four laboratories, 
for an Intercomparison Suite, and three other Lepidolites. (*1.S. was sample a 
part of Intercomparison Suite? ) 


Radiogenic Normal ™S: b Ave Ave 
Sam pli (ppn Sr Sr (x 10°) (10° vr Book 
(ppm) (ppm) units) Diam 
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tories. The sample was very fine grained (0.25mm) and no compositional 
differences between subsamples are indicated by the analyses by the several 
investigators. Three of the four laboratories which analyzed L-N45 reported 
standard deviations for their measurements. These are given in table 3, and 
compared with the actual ranges of the measurements reported, (The ranges 
include the Harwell determinations.) These data indicate that the accuracy 
limits stated by the several investigators are truly indicative of the absolute 
iccuracies of the determinations. 


3 


Relative Standard Deviations (“oS.D.) for Analyses of the Lepidolite 
from Bikita Quarry (L-N45) 


Relative (%) Standard Deviation Actual range 
MILI D.T.M. of values 
Rb Concentratior + 1.6% r12% + 2.8% 4.7% 

Sr Concentration + 2.1% 1.4% 7.9% 


+ 33% +959 + 5.5% 8.5% 


Quantity Measured 


The ages reported for the other Intercomparison Suite samples show 
similarly close agreement. Indeed, in the entire group of 37 isotope dilution 
ige determinations reported for eight localities, only three measurements differ 
from the error-weighted averages of all the measurements made on samples 
from the same locality by more than the standard deviations reported for 
them*, and nearly all of these deviations are less than 6 percent of the values 
reported. 

Concentrations.—Many of the M.I1.T. concentration values reported are 
lower than the corresponding D.T.M. values. We have attempted to determine 
whether these differences are due to differences in spike calibration or pro- 
cedure at the two laboratories, or some other cause. 

In table 4 we have listed the differences between the two sets of analyses, 
with the standard deviations applic able to each value, While there are no sta- 
tistically significant differences in the “Age” column, there are significant 
differences in the Rb and Sr analyses. 

Discussion.—Heinrich, et al. (1953) have reported large variations in Rb 
ontent in lepidolite, Variations of 10 percent are even found from leaf to leaf 
in “book” micas, On the basis of their work, one might have expected concen- 
tration differences in the very coarse Varutrask sample as large as the 13 per- 
cent difference between the M.I.T. and D.T.M. results. But if such differences 
do occur, the perfect agreement between the two fragments analyzed at D.T.M. 
is most remarkable, Still, evidence that variations of nearly comparable mag- 


nitude can occur even in a finer-grained sample is provided by the three 


D.T.M. analyses of the Bob Ingersoll sample: one sampling has about 5-6 per- 


cent less Rb than the others, and comparably smaller **Sr. Hence, it is possible 
that the even larger differences between the M.I.T. and D.T.M. analyses of the 
Pala sample may be real, small-scale sample variations, However, here the 

Leaving out of consideration the very early McMaster determination for L-A106, for 


which (it is now known) too small a Rb sample was taken, so that the reported analysis is 
oo high because of contamination. 
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Taste 4 
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Determine Whether Differences Represent Analytical Errors o1 
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Spattering of solution during digestion could cause such an error, since, 
in our pro¢ edure, although Rb spike was added only to an aliquot consisting 
of 1 to 5 percent of the total solution after digestion, Sr spike was not added 
to the sample in the first instance until it had been digested, We adopted this 
procedure because we were more interested in the age values than in the con- 
centration values per se. But we were aware of the spattering problem and took 
pains to avoid such losses; so, we consider this a very unlikely explanation. 

Nor does it appear that any of the small differences in the procedures 
used at the several laboratories (see section above) could have caused the ob- 
served differences. 

It is our considered opinion that there may be small but real differences 
n spike (or standard solution) calibrations at the several laboratories: pos- 
sibly, analyses of identical samples would run 1-3 percent lower in Rb and 2-4 
percent lower in Sr at M.I.T. than at D.T.M.; also (on the basis of the very 
limited information available) the Harwell Rb calibration appears to agree 
very well with D.T.M.’s, while their Sr calibration appears to be more like 
and several percent different from D.T.M.'s. 

Differences between the results reported larger than those which may be 
due to these possible differences in the calibration standards, and which are 
outside the limits of error stated by the investigators, we at present believe to 
be due to compositional variations within the samples distributed (except that 
the Pala difference may be due to an error in sample identification, as men- 
tioned). However, the weight of the evidence in favor of this interpretation is 
not great, and additional analyses of these samples are greatly to be desired to 
make possible a more certain conclusion. 

Nevertheless, whatever their causes, this study shows that differences are 
not an important source of error in age determination. 

Other samples.—Large differences in Rb and Sr concentration values are 
ilso reported for the pairs of samples collected independently from particular 
sites (table 3), analyzed at D.T.M. and M.I.T.; but here again the radiogenic 
Sr/radioactive Rb ratios and calculated ages for each locality agree with each 


other within the limits of experimental error. We believe that these “sympa- 
thetic” differences in Rb and radioactive Sr concentration also for the most 
part truly record variations in lepidolite composition, However, an alternative 
explanation is the possibility that the samples analyzed contain different 
amounts of impurity minerals of markedly lower Rb content (e.g., calcic 
feldspars, quartz, etc.). But since all samples were routinely purified extensive- 
ly at both laboratories, it does not seem likely that this can be the explanation. 
The possibility that these differences represent analytical errors, is, we believe, 
ruled out by the excellent agreement between the ages determined at the two 
laboratories. 


That agreement in age between independently-collected samples from par- 


ticular localities is excellent in spite of large differences in Rb concentration 
appears also to indicate that lepidolites are relatively unaffected (at least as 
regards Rb and Sr) by surface weathering, The negligible concentrations of 
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nonradiogenic Sr found in these samples show with certainty that no Sr was 


idded to these crystals during their life history. Further, Sr in the slight con- 
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Comparison of isotope dilution, optical spec trographic and neutron acti 
ym analysis.—It is of interest to compare these isotope dilution analyses 
unalyses of the same samples by other methods. 
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5 
Is Normal Sr Found in Intercomparison Suite Lepidolites Mineral Sr 
or Contamination? D.T.M. and M.L.T. Analysis Compared 


Lab Radiogeni Sample Radiogenic Nonradiogenic Calculated 
Sr (ppm) Weight (mg) “Sr JY oft Sr (ygm) normal Sr 
total Sr Absolute concentration 
amount (ppm) 


DTM 356 40,2 

MII 32 1000.0 

DTM 84.7 60.1 

DT 81.6 07.2 

DTM 75.4 50.1 

MIl 80. 1025.0 

DTM 73.8 137 

MII 1000.0 

DTM 56.1 

284.3 0.15 
DTM 103 7.2 0.42 
DTM 30.6 0.26 


MIT 12% 110.5 97 0.43 


his results on a single, synthetic standard which he made up to have an “ace” 
of 1400 X 10° years. (He used a half life of 58 10° years in calculating 


> 


volts while the other (Rb 4202) is an arc 


the Rb/Sr ratio.) As one of the analysis lines used (Sr 4077) is a spark line 
‘ 


with an excitation potential of 8. 

line of 2.0 volt excitation potential, they do not exhibit sympathetic behavior 

during arcing. Indes d. Ahrens was forced to attac h quite large limits of error 

to these analyses, both because of poor sample analysis reproducibility and be- 
use of the possibility of systematic error. 

Dr. Ahrens and his co-workers at M.I.T.’s Cabot Spectrographic Labora- 
tory described improvements in optical spectrographic Rb-Sr age determination 
in several later papers: 

a Ahre ns and Gorfinkle (195lb), re port re pli ate analyses ot a single sample (S.E. 
Manitoba) by an improved method. Sr and Rb are measured separately, each vs, an ap- 
priate internal standard; a half life of 50 x 10° years is used in age computations, 
lhe age found is in good agreement with that reported in reference 2. 
b. Ahrens and MacGregor (1951) report ages for a suite of samples from Africa, in- 
uding several from Bikita District. In this case, analyses are made relative to the 
Manitoba sample. Taking the age of this as 2300 10° years,” absolute ages are calculated 
the African samples 
c. Ahrens and Gorfinkle (195la) describe the potassium internal standard method 
of analysis for Rb used in the most recent of the papers just described. 


To permit comparison of the optical spectrographic ages with the isotope 
dilution results, we have recalculated all of the Cabot Laboratory data using 
the half life value 50 X 10° years. Such recalculation, for example, changes 
the age of Ahrens’ synthetic standard from 1400 X 10° years to 1140 X 10° 
years. The spectrographic ages are plotted versus the 1.D. ages in figure 1, In- 


Which was the average of the then available Rb-Sr and U-Pb data. 
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\ inalyzed by Ahrens, et al. As can be seen, this 
he other ages into agreement with their isotope dilution 
yunterparts h computed averages are plotted as diamonds with standard 
leviation indicated. ) 
lable 6 compares element concentration analyses by spectrography and 
neutron activation, for localities in the Intercomparison Suite, with the isotope 
dilution 
195la, 1951b), using internal standardization techniques; the neutron acti 


vation results are by Webster, et al. (1957). 


lata. The optical spectrographic analyses are by Ahrens and Gorfinkle 


The two neutron activation Rb analyses agree well with isotope dilution 
results tor tl e same sample _ and the single optical sper trographic Si analysis 


is in accept ible agreement with the ¢ omparable isotope dilution analysis, How 
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ever, the spectrographic Rb analyses are all 10 to 20 percent higher than their 
isotope dilution counterparts.'' The difference is greatest for the sample with 
the lowest Rb content. This error is in the direction which will make the optical 
spectrographic “ages” computed using these analyses too small, Significantly, 
optical spectrographic ages reported by Ahrens, et al. (as corrected for half 
life) actually are too small; but, if the S. E. Manitoba (L-A109) value is 
multiplied by the inverse of the ratio of the spectrographic Rb analysis to the 
average of the M.I.T. and D.T.M. isotope dilution Rb analyses,’* one calculates 
an age of 2400 X 10° years, and the “ages” determined for this sample by the 
two methods now agree to within their originally stated errors, Similar correc- 
tion of the other spectrographic ages is not justified because the Rb analyses 
recorded in table 6 were not used to calculate the ages plotted in figure 1, 
which were, as mentioned, based simply on line intensity ratios. 


In spite of the ity of this comparison data, as result of this study the 


writers feel that it is st possible that an optical spectrographic Rb/Sr method 


» dete rmination ¢ ct pt ible accuracy could be developed for lepidolites 
ther minerals containing negligible Sr at crystallization, through the use 


ndards calibrated by sotope dilution and/or neutron activation. Experi- 


LUSIONS 


es of five lepidolite samples ranging in age from 

10° to 2600 years, made in four laboratories (one sample having 
inalyzed by all fi yratories) indicates that reproducibilities of better 
L percent can be obtained routinely for (daughter “Sr/parent *‘Rb) 

» and apparent age. by the mass spectrometric stable isotope dilution tech- 


e determined with a precision better than + 2 


Interlaboratory ! is nearly as good as intralaboratory precision 
lata, indica oncentrations found by isotope dilution analysis 
ite sense, as well as reproducible within a par- 
iracy could be further improved by the ex- 
laboratories for comparison analysis, 
The smallness of the amount of common Sr contamination introduced in 


the processing necessary to this age work (of the order of 0.2 micrograms for 


om OT; mples m ike : possible, he« ause of the high Rb content of lepido- 


lites. the iccurate determination ol the ages of lepidolites as young as l 
10° years. Ancient lepidolites, such as the Bikita Quarry, So, Rhodesia deposit, 


provide copious, inexpensive sources of essentially monoisotopic **Sr, which 


make possible large-scale stable isotope tracer work with Sr at low cost. 
\ comparison of age determinations for particular sites, in which samples 


whose Rb contents differed by as much as 30 percent were used, shows no dif- 


have found errors in the same direction in optical spectrographic determination of 
concentration in two rock standards (9) 

Since one has no clear way to decide which of the I.D. samples of A109 the optical 

spectrographic sample most resembled in composition, the value 0.602% “Rb was arbi 


trarily selected. 


oom 
ue. Both Rb and S ees 
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ferences arising from this cause; hence, there probably has been no leaching 
of parent “‘Rb or daughter “*Sr from these samples, The very low (<1 »gm) 


nonradiogenic Sr content of the samples indicates low original Sr content and 


no addition of Sr during sample history. The apparent ages determined are 
very probably true ages within the stated limits of error of the experiments, if 
the half-life of **Rb is taken correctly at 50 10° years. 
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FOR ALBITE COMPOSITION 
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ABSTRACT The reaction, albite jadeite + quartz, has been studied between 15 and 
25 kilobars, 600 and 1000°C, Two systems of generating pressure were used with mutually 
consistent results: a “simple squeezer’, and a nitrogen system with hydrostatic pressure. 


The equilibrium line is given by the relation, P 6000 (+ 500) + 20 (+ 2) T, 


where P is in bars, T in degrees centigrade. This is in good agreement with the curve of 


Kelley and others based on thermochemical data: jadeite and quartz lie on the high- 
pressure side of the boundary. 

The melting curve of albite has been followed to 15 kilobars; it may be represented 
by the relation, 17 1115 + 0.011P. The rise of melting temperature with pressure is 
smaller than had been estimated from the difference of volume between crystal and glass 
it ordinary temperature 


INTRODUCTION 


Reactions for the formation of jadeite (NaAISi,O,) have been discussed 
by Yoder (1950), Yoder and Weir (1951), Kracek, Neuvonen and Burley 
(1951), Adams (1953), and Kelley, Todd, Orr, King and Bonnickson (1953) 
in terms of the available thermochemical data, which are exceptionally com- 
plete for the minerals involved in these reactions, Retaining the designations 
given by Kracek et al, the reactions in the absence of water are: 

A. Albite (NaAISi.O.) Jadeite (NaAlSi.O,) + Quartz (SiO.) 
B. Albite (NaAISi,O.) + Nepheline (NaAlSiO,) 2 Jadeite (NaAlSi.O, 
Nepheline (NaAISiO,) + Quartz (SiO,) Jadeite (NaAlSi,O,) 

Reactions B and D have been explored experimentally by Robertson, 
Birch and MacDonald (1957). Reaction D, for which the free energy change 
is most favorable to the production of jadeite at ordinary pressure, is not 
realizable because nepheline and quartz first react to form albite, albite plus 
nepheline, or albite plus quartz, depending upon the relative amounts of 


nepheline and quartz, Reaction B was found to have the equilibrium line 


represented by P 1000 + 18.5 T, with P in bars and T in degrees centi- 
grade, between pressures of 12 and 25 kilobars (temperatures of about 600 
and 1200°C). The following report deals with the remaining reaction A, 
albite jadeite + quartz. 

The heats of reaction on which the estimates of free energy differences 
ire based show an uncomfortable amount of scatter, leading to marked un- 
certainty in the expe ted positions of the equilibrium lines. The numbers 
idopted by Adams, Kracek and Kelley for reactions A and B are shown in 
table 1, all for 25°C and 1 atmosphere. The preliminary values of entropy used 
by Kracek and by Adams differ by several hundred cal/mole (of TAS) from 
the final ones used by Kelley et al. The basic data for AH are from Kracek: 
these show differences of the order of 5000 cal, mole, apparently associated 
with the two different albite samples and with the technique of grinding the 
jadeite. Adams takes a mean value for the four combinations given by Kracek; 
Kelley selects the combination: Varutrask albite, Burmese jadeite, For reaction 
A, the individual calculated values of AF (25°C, 1 atm.) range from +70 to 


3320 cal/mole, corresponding to pressures of from 170 to 8000 bars; for 
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trom 780 to 1020 cal/mole, corresponding to +980 to —5S000 

The (extrapolated) intercept at 25° found experimentally for reaction 

was 1500 bars (+1000), which corresponds to AF +1200 (+800) 
cal/mole. This agrees most closely with the values found for the combination: 
Varutrask albite, Japanese jadeite (ground in mullite mortar), though neithe1 
of these materials played a part in the experiments al high pressure, The ex 
perimental determination of the equilibrium, albite jadeite plus quartz, 


ves another opportunity tor correcting the free energy values. 


Tani 
Parameters. at 2 atmosphere 


Material 
Albite Jadeite 


Bu 
70 Jp.ag ile 
1650 Ip.mullite 
860 Ip.(av) 
2350 
1740 
3320 Jp.mul 


fav) 


AH) 


Jp.mullite 
Ip fav) 
Bu 

p.agate 
Jp.mullite 


Ip (av) 


AH) 
Ext apolated 
The I tabil albite, in the temperature-pressure plane. is 


bounded on the high-pressure side by the equilibrium curve for Reaction A, 
on the high-temperature side by the albite melting curve, A few observations 


recorded below determine the position of the melting curve up to about 15 
kilob irs; as the curvature is probably slight, this may be extrapolated without 
serious error to the intersection with the line of Reaction A to define the re 


gion of stability of pure ilbite 


2] 
Uhe 
\ \ 0) 
vr 
i 610 \ Bu 
kx lated 
\ 
102 Jp.agate 
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Va 
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MATERIALS 

The reactants were natural jadeite from Burma (Holden Collection No. 
104644, Harvard Mineralogical Museum), albite from Amelia County, Vir- 
ginia, and quartz from Lisbon, Maryland; the quartz was the purified material 
used by Kracek et al.. secured through the courtesy of the Geophysical Lab- 
oratory. New analyses of the jadeite and albite by J. Ito, Harvard Department 
of Mineralogy, are given in table 2. These are close to the analyses published 
by Kracek et al., showing a small amount of diopside and acmite in the jadeite, 
and a trace of potassium feldspar in the albite. The analyses have been cal- 
culated as water-free, though about one percent of H.O~ was reported in each; 
this is considered to represent water adsorbed on the surfaces of the powdered 
materials. 

Purified Varutrask albite, kindly supplied by J. F. Schairer of the Geo- 
physical Laboratory, was used in a few experiments. Albite glass, prepared by 
Dr. Schairer, was used as a starting material for the determination of the albite 
freezing curve. 

2 
Chemical Analyses of Reactants 
(calculated as water-free) 


3 1 5 6 
SiO 9.3 ; 59.45 68.18 68.57 68.82 
rio 
ALO i : 25.22 19.83 19.86 19.51 
8 35 16 03 13 
FeO ) 02 
MeO 8 04 
14.38 11.65 11.19 
2 ke .29 
100.00 
l Jadeite, Burma. Analyst J. Ito 
. Jadeite, Burma. Analyst E. G. Zies 
3. Jadeite, theoretical, NaAlSisOc. 
Albite, Amel a, Virginia. Analyst Ito 
». Albite, Amelia, Virginia, Analyst E. Chadbour: 
6. Albite, Varutrask, Sweden, Analyst E. G, Zies 
Albite, theoretical, NaAlSi,O. 
EXPERIMENTAL DETAILS 

The procedure follows closely that described by Robertson, et al, and by 
Clark, et al. (1957) with a few later improvements. A water jacket around the 
supporting rings of the pressure system keeps the temperature at the pressure 
gauge below 50°C, eliminating a troublesome correction in the pressure meas- 
urement. A Minneapolis-Honeywell Pulse-type Pyr-o-vane controller has been 
used to control the furnace temperature; a shunt around this controller carries 
most of the heating current, the controller affecting only enough to prevent 
long-time drift of temperature. Reduction of length of the platinum capsules 
containing the charges from 58 to 44 inch has reduced the difference of tem- 
perature between the ends of the capsules. The smaller size of plotted rectangles 
indicating variation of temperature and pressure during the exposure of 
samples shows the effectiveness of these arrangements. 


Except for a few runs in which silicic acid was introduced in place of 
quartz, the ingredients were dried by baking in the capsules at 300°C for 24 
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hours They were then seale d by welding. Silic ic ac id had little effect on this 


reaction in the hydrostatic system but appeared to promote reaction in the 


squeezer, In general the crystals produced were exceedingly small. 
| | 


The equilibrium curve was sought from both sides, by causing albite to 
break down to jadeite + quartz, and by causing charges of jadeite + quartz 
to form albite. It is, of course, impossible to reach a particular combination ol 

h pressure and high temperature without passing through intermediate 
values of pressure and temperature, and the possibility that a reaction may take 
place at some intermediate and unknown point has to be guarded against. The 
ipproach to the final values is made as rapidly as possible, and it is made 
through the field of stability of the starting material. When starting with albite, 
for example, the temperature is raised first, then the pressure; when starting 
vith jadeite plus quartz, the pressure is raised first, then the temperature, I! 
no change is observed, no importance ts atta hed to the result: only changes ol 
phase have been tabulated and plotted. These occurred in about one-third of 
the runs 

Identification has depended upon the Norelco X-ray diffractometer, as 

1ided by the microscope whenever the end products were sufficiently 
A small amount of albite is readily detected against a background ol 
jadeite and quartz lines by the presence of the strong lines at about 24 


JADEITE + 
QUARTZ 


ALBITE 


@ JADEITE + QUARTZ 


PRESSURE 20 KILOBARS 30 


tem perature pressure plane. Rectangles designate 
em, with dimensions corresponding to variations of 
Circles (arbitrary size) designate determinations in 


bd 
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(Copper radiation). The detection of a small quantity of jadeite against 
a strong albite background is much more difficult, as there are albite lines 


close to most of the normally strong jadeite lines; extension of the diffractome- 
ter pattern to include the jadeite line at 26 13.7° is helpful, no albite line 
interfering at this point. This is only the third or fourth strongest jadeite line, 


however, so that a small amount of jadeite may still not be recognized, The 
pressures at which albite is found to break down to jadeite and quartz may 
thus be slightly higher than the minimum pressures required; on the other 
hand, the reaction of jadeite plus quartz to form albite, relatively easy to de- 
tect, sets a lower limit for this pressure which is not far below the observed 
points of albite breakdown, There is a margin of uncertainty of about 1 kilo- 
bar which probably results from lack of precise control of the pressure- tem- 
perature conditions. 

The equilibrium could not be followed below about 700°C with the hydro- 
static system. Some additional experiments were made at 600°, 700°, and 
800° with a “squeezer” of the type described by Griggs and Kennedy (1956). 
As shown in figure 1, the consistency between the hydrostatic runs and the 
squeezer runs is good. The range of pressure within the squeezer charge is 
probably of the order of several kilobars unde sood conditions, but possibly 
much greater on occasion, With average pressures spaced at intervals of 5 
kilobars, the results show little confusion for the present system. 

For the investigation of the melting curve, two capsules, one containing 
ulbite glass, the other crystalline albite, were heated under pressure to a final 
temperature, held here for some hours, and quenched. It had then to be de- 


3 
Experimental Results, Hydrostatic System 
albite: Jd jadeite; Qu quartz: Si silicic acid: gl class ) 


Time Pressure Temperature 
hours bars ( arg Products 


19.850-20,040 993-1003 

20,440-20,500 1012-1023 

19.760-19,840 1069-1080 

15.720-15.800 848-866 

15.000-15,000 846-854 

15,080-15.200 840-856 

15,160-15,120 800-832 

16900-17000 796-803 

16.930-17 000 744-755 

2?,720-22.840 738-761 

21,280-21,760 840-863 

23,400-23 440 799-807 + Qu 

2? 400-22 640 798-807 +Qu+ Ab 
+ Qu 

1). 000-20,120 794-810 me +Jd+Qu 


+ Qu 


Jd+Qu+ Ab 
23.480-23.520 921-939 + € +Jd+Qu 
+ Qu 


24,320-24,730 890-904 


760-24,880 994-1020 


+ Qu 


Ab 
) 6.0 
6 
87 
15 
LOO 
103 
0 14.0 
ld+Qu-+ Si 
109 4.5 CSC Ab Ab 
125 11.75 Id+Qu Ab 
Ou+ Ab(5%) Ab 
29 0 
Ab-+S \b 
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TABLE 4 


Experimental Results, Squeezer 


Ab albite; Jd jadeite; Qu quartz; Si silicic acid: gl 


rime 
hours D arg Products 


20 10.000 7 Ab+Qu 
1.0 12,000 7 t Jd+Qu- 
? 15.000 8 Ab-+ Jd t 
10 10.000 8 +Jd- 
? 10,000 8 + Qu- + Qu 
L5 15.000 } Qu 
0.5 8,000 8 + Id- 
1.0 5.000 7 + Jd- 
1.0 20.000 7 + 
1.0 5.000 70) d+Qu-+ 
1.0 0,000 8 1+ +Jd-+ 
18.000 } Id T 
OO) +Qu-+ + Qu-+ 
5.000 ) t Qu t 
) } d+Qu-4 
+ Qu- 
+ 

+ 
+Qu-+ 

Qu 
Ou 


5 
ts Near Melting Curve 
ulbite: el class) 


Dime Temperature 


hours 


1198-1210 
1240-1260 
1274-1279 
1298-1314 
1135-1154 


1.20 000 1341-1366 
UK) 1386-1395 
880-22.960 1321-1358 
9 800-20.040 1286-1323 


ed in proof; 


rask albite. 
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Ab(?) 
Ou 
Qu 
Qu 
Qu 
Ab( 
Ab 
) Ou 
Qu 
| Ab( 7) 
14 Ab 
Abi 
4 Ab 
Ab 
| Qu 
64 Ab 
Ab 
Ou 
\| 
\ 
(ha Products 
Ah ol 
4.0 Ab Ab 
Ab Ab 
Ab 
Ab ¢ gl 
zl Ab gl 
bs Ab gl Ab 
Abs gl 
) Ab gl 
Ab gl gl 
Ab gl+-(Ab 
0.2 Abt (Ab 
134 4 Ab? Ab 
l 2.0 Abi Ab 
Varu 
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termined whether the glass contained crystals, or the originally crystalline ma- 
terial contained glass. The properties of albite glass are relatively suitable for 
this procedure, since it crystallizes readily enough under pressure, but not so 


readily that it cannot be quenched as a glass. 


RESULTS 

The individual experiments are summarized in tables 3, 4, and 5, and 
plotted in figure 1, Squeezer runs are designated by circles of arbitrary di- 
ameter, hydrostatic runs by rectangles whose dimensions correspond to the 
range of temperature and pressure observed during the high-temperature ex- 
posure. Only experiments where the charges definitely changed have been 
retained, though there are many where negative confirmation was obtained. 
Thus a number of runs in the albite field where albite was unchanged have 
been omitted, and similarly for the jadeite + quartz field. 

The experimental equilibrium line is drawn with a slope of 50°/kilobar 
or 20 bars/ degree. This is derived from dT/dP = AV/AS with an uncertainty 
of perhaps 10 percent.’ With the relatively wide margin between the experi- 
mental points in the albite field and those in the jadeite + quartz field, the 
uncertainty of the slope based on these points is perhaps 20 percent. If the 
slope is assumed to be fixed by the thermochemical data at 20 bars/degree, the 
intercept at O°C cannot differ by more than about 500 bars from 6000 bars. As 
equation of the equilibrium line, we take P = 6000 (+ 500) + 20 (+ 2) 
lr, with P in bars and T in degrees centigrade. 

Che equilibrium line calculated by Kelley, et al.. is shown in figure 1 as a 
broken line: this is in good agreement with the experiments, showing no sig- 
nificant difference if allowance is made for the uncertainty of the thermo- 
chemical data. This agreement suggests that the rather large difference be- 
tween experimental and calculated curves for reaction B (albite + nepheline 

2 jadeite), if we use the average AH. might arise from the heat of reaction 
Oi nepheline, rather than from that of albite o1 jadeit . as assumed by Robert- 
son, et al, The experimental lines can be extrapolated to 25°C to give values 
of AF, and thence of AH, with the aid of Kelley’s values of AS; these figures 
are included in table 1. It seems probable that this method of finding AH for 
these reactions is more reliable than the method of taking the small difference 
between relatively large heats of solution. 

The calculation of the rise of melting point with pressure based on the 
difference of volume between crystal and glass at ordinary temperatures 
(Goranson, 1938, p. 82) gave 26°/kilobar. The observed change, fairly well 
defined by the measurements, is 11 + 2°/kilobar. Similar discrepancies have 
been noted for other melting curves '( Yoder, 1952; Clark, 1958). Yoder’s sug- 
gestion that there is an abnormally large thermal expansion of the solid just 
below the melting point accounts very well for the usual sign of the discrepan- 
cy, though there seem to be no direct observations of this: effect. The new 
value for albite is close to, but slightly lower than, Yoder’s value for diopside, 
13°/kilobar to 10 kilobars. (See also Boyd, 1958, for the melting of diopside 


it higher pressures.) 


With AV 17.5 cm*/mole and AS 8.3 cal/deg.mole, dT/dP 50.4 deg/kilobar, 
i] 


IP 19.8 bar/deg 
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These changes lead to modifications of the tentative diagrams given by 
Robertson. et al. for the relations in the system NaAlSiO,-NaAlSi.Oxg, assumed 
to be binary, at 15 and 30 kilobars. The albite melting point must be reduced 
and the position of the lines defining the reaction, albite jadeite + quartz, 
corrected. Some further investigations of melting of albite-nepheline mixtures 
indicate that the nepheline + liquid field found by Greig and Barth (1938) 
for mixtures of jadeite composition retreats, as the pressure is increased, to- 
ward the nepheline side, so that at 15 kilobars, a nepheline-albite mixture of 


jadeite composition is in equilibrium with albite plus liquid, and at a tem- 


perature somewhat higher than is shown on the diagram, Correction of these 


diagrams requires more observations than are presently available, and is left 


for a later occasion 


NTINENTAL” 


reactions plotted as function of d pth, with 


Figure 9 of Robertson, et al. has been revised and appears as figure 2 of 
the present paper. The breakdown of pure albite occurs at pressures greater 
than those of the base of the normal continental crust for temperatures greater 
than 200°C, The difficulty noted in the earlier paper with the interpretation of 

sub-oceanic discontinuity as a phase change involving jadeite of course 

sists. The experimental results are still too incomplete to afford a satisfactory 


basis for discussion of the reaction. gabbro eclogite. 


4 

3 
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ABSTRACT An s made of the various possible methods of investigation in the 
Earth sciences, It is shown that, whereas much effort has been directed toward the obse: 
, the theoretical analysis of these data and that of the intuitive inferences 
1 therefrom has be hat lacking, Hope is expressed that mathematical analys 
Ka sciences as it has in othe: discipline s. 


After publication of the writer's book on geodynamics (Scheidegger, 


comments have appeared in the literature 


(MacDonald, 1959; Beloussov, 1959). It appeared that these comments were 
not concerned wil ‘cific faults in the book, but rather that they implied a 
he writer's method of investigation, This method is on 

irguments originating in mathematical physics to 


the Earth sciences, Although such a procedure is by no means new (the great 
classical works of Darwin, 


of mathematical analysis in the Earth sciences is still the subject of hot con- 
troversy. 


This cont 


seemed to him, 


1958). several unfavorable 


wholesaie rejects 


of systematical 


Love and Jeffreys bear witness to this), the validity 


th 


he writer partly due to misunderstanding. It 
1 somewhat closer analysis of the various meth- 
the Earth sciences and of their relationship with one 
another might be in order. It was hoped that this would help to untangle some 


1as arisen and, possibly, to undo some of the discredit 


il method has fallen. 


h this plan, let us first review the various methods 


ial 
ods of investig 
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~ disposal, Any science starts with ob 
1. This is a basic tenet whether the science is 
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trical field or with the structure 


nm ¢ 


simply collec ted and des ribed. 
1 can be accumulated, Such informa 
$s, Maps, or ve rbal des ription. The body of 


s the first and basic stage of scientific research 


ture 


the information, it is necessary to bring som: 
y class 


ifying the various observations into cate 
ouped into families and genera, sedimentary strata 
iples of stratigraphy, rocks are ordered a 


ul composition, In this fashion, one ends up with a 


Lo he prin 


stem il lected information becomes easily accessible and new 
information car issified into its proper category. Thus, a basic 
ted which must remain the foundation for any 


body ol knowle dae 


lurther progress 


umulat 


[he attempt will often be made to set the qualitative descriptions referred 
to above on a more quantitative basis. Many properties of the features investi- 
in tact, be 


e measured, In the Earth sciences, it has most- 


I 
cated qualitatively 
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ly, but not exclusively, been the geophysicists who have made use of measuring 
devices. The gravity field, the depth of the Mohorovicié Discontinuity, terres- 
trial heat flow and many other quantities can readily be measured, Thus, to 
the earlier, purely descriptive information there is now added a set of data 
that can be expressed in terms of numbers. 

After a body of information has thus been accumulated, the scientist will 
want to draw inferences. First of all, this can be done in a rather direct way. 
Taking a particular body of data (which may stand apart from other data), 
one can ask oneself how this body of data is explainable. In most instances, 
this requires constructing by intuition a hypothesis of a process that might 
have produced such data. Taking an example from the Earth sciences, one may 
look at paleomagnetic data and one might infer that continents drifted over 
large distances during geologic history. Similarly, if one were to consider the 
nappe structure of the Alps, one might infer that crustal shortening is playing 
a rather important rdle in Earth processes, Unfortunately, the inferences or 
conclusions thus drawn are always tied up with a hypothetical process and the 
processes inferred from different sets of data are often mutually contradictory. 
Thus, the question arises as to how the different hypotheses that intuition has 
invented can be reconciled, A natural way to go about this is to try to collect 
more evidence which would have a bearing upon these processes. Unfortunate- 
ly, it lies in the nature of the Earth sciences that no direct observations of these 
processes can usually be made; the latter may take many millions of years for 
their completion, which is much too long for any human being to witness. In- 
direct methods are therefore required. 

It is here that a systematic analysis of the physics and chemistry of the 
processes involved can provide a valuable tool. This method has often been 
termed “mathematical geology”, because physics and chemistry, being “exact” 
sciences, use mathematics as the chief means of expressing their findings. 

Let us now have a more specific look at the “mathematical” method in 
geology, i.e. at the method of testing the various geological hypotheses with 
regard to their physical and chemical contents. In applying, say, physics to 
eeology, the use of mathematics is mandatory, Physics is a science whose find- 
ings can be most easily expressed by mathematical equations. This, at once, 
poses a problem in connection with its application to the Earth sciences: 
Geology, being essentially a descriptive science, does not easily lend itself to 
an expression in terms of numbers. Numbers, however, are the basis of equa- 
tions. 

Thus, in applying physics to geology, it is first of all necessary to express 
geological findings in terms of numbers. This will not often be possible without 
a certain amount of arbitrary systematizing. Thus, it is quite obvious to any 
one looking at a map that the strings of islands off the northeast coast of Asia 
form rather regular chains. One would expect that there ought to be a reason 
for their existence. Many intuitive and mutually exclusive hypotheses for their 
existence have indeed been put forward, and in order to choose amongst them, 
a test of them in the light of physics might prove helpful. In order to do this, 
one must first express the shape of the arcs in terms of numbers, One might, 


for instance, regard the arcs as circles and systematize them into a pattern of 


with which to do geometrical exercises. This was done, for instance, 
by Wilson (1952). Critics of such a systematization usually maintain that a 
procedure of this sort represents too mut h of a tour de Jorce and that one 
should, therefore, confine oneself to a mere description of various disjointed 
facts. However, in this manner it will never be possible to proceed to any 
reasonable test of the various intuitive hypotheses that have been put forward 
istence of these facts. The reason that circles may conveniently be 
description of island arcs is merely that circles are simple geo- 

ures resemblir somewhat the observed shapes. If one can rational 

xplain the circles, features that resemble circles can then be explained by 

ynditions that resemble those that cause circles. 

Since the systematization of data is to some extent arbitrary, different 
systematizations will be possible on the same set of data, Thus, the 

of islands might | garded as parts of X’s rather than as circles 
Wilson’s (1952) interpretation in terms of circles is still the most 

plete that is available to date. Similar remarks could be made with regard 

to many ‘r attempts to dese ribe geology in terms of numbers. If one con- 

demns one type of systematization, the problem is to find a better one, 

fter havin tablished a set of numbers as characteristic of certain 

i features, it is then possible to bring physics into the picture, For a 

sts have been hypothesizing about the origin of the 

\ multitude of “theories” has been invented that were 

These theories contradict one another and 

les of faith rather than of reason, with the proponents 

y a degree of intransigence not usually found in the 

mnfusion would appear to be a careful and syst 

hypotheses in the light of a discipline,-physics, 

tily applied to these hypotheses and theories 

of the theories out as being unreasonable 


»ward that which is “true” by a process of elimi 


It must be understood, thus, that the mathematico-physical method will 


to ‘solve’ problem ol ceology Therefore. if Beloussoy states 


ol 


whether or not the poles have migrated is 

storical geological data, but on mathematical 
on of the viscosity of the Earth's crust. 


bears witness to a common type of misunderstanding of the mathematical 
method, The argument should run as follows: Certain geological and geo 
physical observations suggest that the poles have wandered, If this interpreta- 
tion of observational data is correct, an upper limit for the effective viscosity 
of the Ear ve deduced by the mathematical method, This value, then. 
must be consistent with deductions of the same quantity from other sources of 
information, If it is not consistent, this does not necessarily mean that polat 


wandering did not take place, as it could also signify that the previous inter- 


pretation of observational data leading to the adopted viscosity value was 


iulty, 7 hus. Oo! ly the consistency ol inferences from different sets of data can 
be checked by the “mathematical method”. In particularly clear-cut cases, an 


rrogenetic theory can also be shown to be inconsistent within itself without 
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further reference to other data. This, it appears to the writer, would be a good 
thing to know. 

Many Earth scientists have restricted themselves to carefully collecting 
evidence about the Earth, an endeavor which is basic to all further studies. 
\lany have proceeded to draw all kinds of inferences which such evidence 
would allow them to draw, However, it appears to the writer as strange indeed 
that (as Beloussov and MacDonald seem to imply) the only check on such in- 
ferences should come from collecting more evidence, and that it would not 
even be permissible to investigate the self-consistency of such inferences in the 
light of related sciences, such as physics and chemistry. Surely, if we are to 
remove the speculations about geodynamics from the realm of mysticism and 
dogma into that of reason and enlightment. no means. not even mathematical 
physics, ought to remain untried. 

Phe mathematical method in geology, thus. would appear to be able to add 
1 valuable store of knowledge to the Earth sciences. It is unfortunate indeed 
that its function is so often misunderstood, Naturally, up to the present time, 
only a mere beginning has been made in investigating the mathematico- 
physical significance of many geological observations. The mathematical meth- 
od has fallen into discredit by many gross errors that have been committed at 
times. Nevertheless, in spite of all these short-comings, the writer would like 
to hope that the understanding of the power of mathematical analysis will grow 


in the Earth sciences as it has in other disciplines. 
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REVIEWS 

Kritik und Fortbildung der Relativitatstheorie; edited by KARL SAPPER. 
P. 281. Graz, 1959 (Akademische Druck- Verlagsanstalt, $11.50) .—This 
book is a miscellany of articles by authors who are philosophers, or physicists, 
or both, Seven of the papers are in German, four in English, Their purpose is 
to examine the philosophical meaning of the special theory of relativity. 

The gist of the book will become clear from a discussion of four contribu- 
tions which are selected largely at random. Professor A. Griinbaum discusses 
the philosophical foundations of the theory and presents what may be termed 
an explanation of the ideas proposed by Reichenbach in his well-known books”*. 
Griinbaum’s paper contains many interesting insights and establishes connec- 
tions with all the m 1jor interpretations of relativity which have been proposed. 

The article by P. Moon and D. E. Spencer is more analytic in nature and 
tries to establish relativity on a new physic al basis. These authors affirm that 
there are two ways of accounting for the Michelson-Morley experiment, The 
first is to keep Maxwell's equations intact, but to reject the Galilean trans- 
formation in favor of th« f Lorentz; the second is to keep the Galilean 
transformation but modify Maxwell's equations, Einstein chose the former 
course, Moon and Spencer propose the latter. The authors seem to be unaware, 
however, that their approach, which is known to most physicists, destroys en- 
tirely the simplicity and the conceptual beauty of the theory of relativity, 

Dr Sapper examines \ ith great care the various meanings of the word 
‘relativity’. After making some interesting distinctions, he concludes in some- 
what heterodox fashion. that the constant C, usually regarded as velocity of 
light in vacuo, is conceptually different and requires an interpretation in terms 
of measurement. The science of relativity, he holds, lies in the relativity of 
measured values, not in any properties of the objective process of nature 
themselves. 

A. Wenzl, the distinguished philosopher of the University of Munich, 
makes a contribution which is particularly noteworthy when viewed against 
the background of the current controversy on the meaning of quantum me- 
chanics. lie resolves what he regards as the paradox of simultaneity (i.e., the 
lack of difference in different reference frames) by suggesting that the in- 
fluence spreads through space in a manner subject to the Lorentz transforma- 
tion is not the physical agency called light, but a “potency”, a probability 
which can be actualized, “corpuscularized”, as the author says, when a meas- 
urement is performed. By this device he wishes to save pre-relativistic intuition 
in connection with actual occurrences in time and space. 

The book is full of novel ind unorthodox points of view and makes re- 
freshing reading HENRY MARGENAU 


(uaternary Stratigraphy and Climate in the Near East: by K. W. Burzer. 
P. 157; 16 figs., 4 maps. Bonn, 1958 (Bonner Geographische Abhandlungen, 
Heft 24, 12 DM) .—The study of Pleistocene events beyond the glacial border 
encompasses many segments of natural science. Geomorphology, paleoclimatol- 


ogy, physical geography, paleontology, archaeology, and biogeography are 


* Reichenbach, X., 1924, Axion ler Relativistischen Raum-Zeit-Lehre: Braunschweig, 
Vieweg; 1928. Philosophie der Raum-Zeit-Lehre: Berlin, De Gruyter. 
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among the sources Butzer draws upon. He is particularly concerned with the 
first two fields. In part his conclusions are summarized in a map of Wiirm 
vegetation zones. Temperate mixed and coniferous forests cover Anatolia, 
mountain tundra occupies the highest plateaus and peaks, subtropical forest 
and scrub borders the western Mediterranean coast, gallery forest follows the 
Nile, and in Lybia the desert endures, despite marginal invasion by steppe and 
semidesert. 

Butzer is bold enough to attempt such a map in the absence of a rich 
paleobotanical record, On the basis of 600-800 meters of snowline depression 
in Iran and Anatolia he infers a maximum temperature depression of 4°C in 
the Middle Wiirm. On the basis of stalagmite development in Mediterranean 
caves he appraises the Early Wiirm as the main period of increased rainfall. 
Apparently this was insufficient to generate true pluvial lakes in Libya and 
Egypt. Finding no fossil pluvial drainage systems in the greater part of the 
Libyan Desert, he assumes that it was not obliterated by more mesic vegetation 
zones. He regards the pluvial periods as relatively brief and temporary inter- 
ruptions of the normal climate of this region. 

American stratigraphers, accustomed to seeking a Xerothermic period in 
the postglacial sequence of the arid Southwest may be surprised to note that 
the period from 7000 to 4400 BP is considered a subpluvial. Terra rossa and 
stalagmites formed in Hagfet et Tera Cave, Cyrenaica; rock pictographs fea- 
ture a savanna fauna in parts of the Libyan Desert; arroyo aggradation oc- 
curred at Wadi Diglah. Butzer notes a parallelism between moisture trends in 
Europe and in the Near East; in Atlantic time both areas appear warmer and 
wetter than at present, in the Sub-boreal both were drier. 

Discounting the theory of Dubief regarding a shifting jet stream, Butzer 
relates both the Near East pluvial and postglacial subpluvial periods to intensi- 
fied Mediterranean westerly circulation. In fairness to Dubief, it may be noted 
that both Ptolemy’s weather record at Alexandria and the relict colonies of 
tropical reptiles in Egypt and Barbary (Braestrop 1947, Vidensk Medd, fra 
Dansk. Naturh. Foren., Bd, 110) seem to evidence past periods of summer 
monsoon rainfall north of the Sahara. Butzer’s case rests partly on diminishing 
evidence of subpluvial conditions southwards from 30° to 23° N, lat. 

While Butzer omits Shalem’s notable dissent regarding the importance of 
Pleistocene climatic change in the Middle East, | assume that his 336 ref- 
erences cover most of the pertinent literature on this region. Since this is the 
first English volume in the Bonn series, some alarming syntax and spelling 
may be forgiven. Readers with a general interest in prehistory and climatic 


change, as well as Pleistocene specialists, will find Butzer’s account highly in- 
formative and stimulating. PAUL S, MARTIN 


Coastal and Submarine Morphology; by ANpre GUILCHER (translated by 
B. W. Sparks and Rev. H. W. Kneese). P. 274; 40 figs., 8 plates. New York, 
1958 (John Wiley and Sons, Inc.), $6.50).—This is a new edition as well as 
a translation of a book first published in France in 1954 under the title “Mor- 
phologie Littérale et Sous-marine.” Being only 274 pages in length it is not an 
exhaustive survey of coastal and submarine geomorphology, however a re- 
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markable amount of information is contained within its covers. It should be 
f considerable value to the American student of this field because most of the 
examples of coastal forms are taken from Europe and Africa. The bibliogra 
phies given at the end of each chapter, include some 800 titles and appear to 
represent a very thorough survey of European literature, although the Ameri 
in publications in the field are not slighted. A number of the references are 
innotated 
It is interesting t O that in his discussion of shoreline movements 
Guilcher iccepts Pleistocene eustatic terraces at 5. 15 to 20. 30-35. 55-60 and 
80 to 100 meters. He believes that the Pleistocene is too recent and short to 
produce a significant number of terraces caused by diastrophic movement. He 


follows Shepard in classifying coasts as primary or initial and secondary or 
sequential, Primary coasts are ria, fjord, glacial lowland, unglaciated lowland, 
ind coasts dominates 
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